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THREE-DIMENSIONAL TRANSONIC FLOW THEORY 
APPLIED TO SLENDER WINGS AND BODIES 
By Max. A. Heaslet and John R. Spreiter 


SUMMARY 


The present paper reexamines the derivation of the integral equations 
for transonic flow around slender wings and bodies of revolution, giving 
special attention to conditions resulting from, the presence of shock waves 
and to the reduction of the relations to the special forms necessary for 
the discussion of sonic flow, that is, flow at free-stream Mach number 1. 
In the vicinity of the body, the disturbance field is then shown to con- 
sist of a two-dimensional disturbance field exten ding laterally and a 
longitudinal field that depends on the streamwise growth of cross-section 
area. This result extends Oswatitsch's equivalence rule to lifting cases, 
provided the angle of attack is small relative to the thickness ratio. 

The correctness of the analysis is checked by examination of Yoshihara 1 s 
numerical solution for sonic flow around a slender, circular cone-cylinder 
and this solution is checked, in turn, by comparison with experimental 
results of Solomon. An example is presented in which the general result 
is applied to calculate pressure and integrated forces on a family of 
slender, elliptic cone-cylinders. An expression is derived which permits 
the ready calculation of the difference in drag of two slender bodies 
having the same longitudinal distribution of cross-section area. Classes 
of wings and bodies are described for which the difference in drag is 
zero and the Whitcomb area rule applies. Experimental data for such a 
family of wings of rectangular plan form are examined and it is shown 
that theory and experiment are in good accord. 


INTRODUCTION 


The equations governing transonic flows are known and well estab- 
lished by favorable comparisons with experiment (see ref. 1 for a resume). 
The difficulties arising as a result of the nonlinearity and mixed char- 
acter of the differential equation for the potential, however, have pre- 
vented the rapid advancement of the analysis such as has occurred in 
recent years with both subsonic and supersonic theory. This is particu- 
larly true for three-dimensional transonic flows and, as a result, perhaps 
greater than usual effort has gone into the determination and utilization 
of relations between solutions. The first of these to be advanced was 
the transonic s imilar ity rule which pertains to the pressures and forces 
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on affinely related wings (refs- 2, 3, and and bodies of revolution 
(ref . 5) • A second relation is the area rule established empirically by 
Whitcomb (ref . 6) which states that ’’near the speed of sound, the zero- 
lift drag rise of thin low-aspect-ratio wing-body combinations is pri- 
marily dependent on the axial distribution of cross-sectional area normal 
to the air stream." A third relation is the equivalence rule of Oswatitsch 
(refs. 7 and 8) which may be stated as follows: "The solution for tran- 

sonic flow around a thin, nonlifting, low-aspect-ratio wing can be obtained 
from that for a nonlifting body of revolution having the same longitudinal 
distribution of cross-sectional area by superposing the difference between 
the two-dimensional harmonic crossflow solutions for the two bodies." 

The area rule and the equivalence rule are, obviously, closely related. 
Further effort needs to be expended, however, in establishing the gener- 
ality and range of validity of these relations and in exploiting the 
results in specific applications. The present paper is concerned with 
this task. 

The problem will be approached through application of the classical 
method of singularities. This is one of the oldest and most fruitful 
methods for solving partial differential equations and has reached a high 
state of development in lin earized compressible-flow theory. There is 
also a considerable body of literature in which the method is applied to 
nonlinear compressible-flow problems by considering the solution of the 
linearized equations to be a first approximation, and iterating to obtain 
second and higher order approximations. The results so calculated are 
good approximations to pure subsonic flows or to pure supersonic flows, 
but it is now generally agreed that the series representation of the solu- 
tion does not converge in the transonic range. Approximate calculations 
by Oswatitsch (refs. 9 &ud 10) have indicated the possibility, however, 
that the method of singularities might be applied successfully in the 
transonic range if the idea is relinquished that the linear solution is 
necessarily the first approximation in the transonic range. This idea 
has been pursued further in references 11 through 15 in which a number of 
improvements are introduced and the results of numerous specific calcula- 
tions are shown. Although the basic equations are derived for three- 
dimensional flow in the latter references, all applications axe to two- 
dimensional flows. The values of the free-stream Mach number, moreover, 
are restricted to values no greater than unity. 

The same general approach has been applied to three-dimensional tran- 
sonic flow around slender wings and bodies by Oswatitsch and Keune (ref. 8) 
and by Harder and Klunker (ref. 16) . In these applications, the principal 
aim is not to determine actual solutions but to derive relations between 
solutions for various bodies having the same longitudinal distribution of 
cross-section area. These two analyses axe not entirely satisfactory in 
a number of particulars , not the least of which is the omission of all 
considerations of shock waves in the body of the analysis. A more impor- 
tant, although perhaps more subtle, point concerns the treatment of the 
cumulative effect of the nonlinear term of the transonic equation when 
the free-stream Mach number Mq is -unity. Harder and Klunker argue that 
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the effect on the induced flow field is small because the term itself is 
everywhere small. Actually, however, the cumulative effect of this term 
leads to infinite contributions at Mq = 1. Oswatitsch and Keune consider 
the cumulative effect but circumvent the difficult by introducing rather 
arbitrarily selected values for Mach number so chosen that the value of 
unity is never inserted into the vital integrals. It is the initial con- 
cern of the present analysis, therefore, to re-examine the derivation of 
the integral equations for transonic flow around slender wings and bodies 
of revolution, giving special attention to conditions resulting from the 
presence of shock waves and to the reduction of the relations to the 
special forms necessary for sonic flow. In contrast to references 8 
through 15, which are concerned exclusively with cases in which the free- 
stream Mach number is no greater than unity, equations are also derived 
herein for the case where the free stream is supersonic. These equations 
are likewise reduced to the special form associated with sonic flow and 
the results are shown to be identical to those which arise from a consid- 
eration of flows with a subsonic free-stream velocity. 

Following the establishment of the basic integral relations for 
transonic flow, special attention is directed toward the case where the 
free-stream Mach number is unity. Here, the integral relations are 
simpler in character, although still nonlinear. Application of a con- 
vergent iteration process leads to the conclusion that the solution for 
the potential has a particularly simple form in the vicinity of the body; 
in common with linearized slender-body theory, the disturbance field 
consists, to a given order of error, of a two-dimensional disturbance 
field extending laterally and a longitudinal field that depends on the 
streamwise growth of cross-sectional area. This result extends 
Oswatitsch' s equivalence rule to lifting cases, provided the angle of 
attack is small relative to the thickness ratio. The correctness of 
the analysis is checked by examination of Yoshihara's numerical solution 
for sonic flow around a slender, circular, cone-cylinder given in refer- 
ence 17, and this solution is checked, in turn, by comparison with experi- 
mental results of Solomon given in reference 18. The results yield a 
simple means of determining the pressure distribution on an entire family 
of slender wings and bodies having the same longitudinal distribution of 
cross-sectional area when the pressure distribution is known for any 
member of the family. Starting with the known solution for the circular 
cone-cylinder, an example is presented in which the general result is 
applied to a family of slender elliptic cone -cylinders. This example, 
which was discussed briefly in reference 1, is examined in detail. It 
is shown that the lift and the load distribution are the same as given 
by linear theory, confirming the ideas advanced in reference 19. Con- 
trary to Whitcomb's area rule, however, the drag depends significantly 
on the cross-section shape. Both the drag and lift of the thin elliptic 
cone -cylinder are shown to be in accord with the transonic similarity 
rules. A momentum analysis of the sonic drag of slender bodies in gen- 
eral is then undertaken and an expression is derived which permits the 
ready calculation of the difference in drag of two slender bodies having 
the same longitudinal distribution of cross-section area. This result 
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confirms the drag variation calculated for the elliptic cone -cylinders 
by integration of the surface pressures. Several large and significant 
classes of •wings and bodies are described for -which the difference in 
drag is zero and the Whitcomb area rule applies without modification. 

One of these is a family of affinely related wings. Exper ime ntal data 
from reference 20 for such a family of wings of rectangular plan form are 
examined and it is shown that theory and experiment are in good accord, 
provided the product of aspect ratio and, cube root of the thickness ratio 
is, in this instance, less than about unity. 

The final section of the report has been written so as to be as 
self-contained as possible and readers concerned solely with applications 
of the theory may find this section sufficient for their purposes. The 
initial sections of the report have been written for readers concerned 
with a more complete understanding of the derivation and limitations of 
the general theory together with the evaluation of the order of error 
incurred in the slender body approximation. 


FUNDAMENTAL EQUATIONS AND BOUNDARY CONDITIONS 


The basic equations necessary for the discussion of inviscid tran- 
sonic flow consist of a set of partial differentia], equations relating 
the velocity components and their gradients at every point, together with 
an auxiliary relation for the velocity jump through a shock wave. For 
thin wings and slender bodies inclined at zero or small angles of attack, 
the differential equations can be simplified if the shock waves are 
assumed sufficiently weak that the flow is irrotational and isentropic, 
thereby permitting the introduction of a velocity potential The 

further assumption of small disturbances leads to the use of a perturba- 
tion velocity potential q> which in Cartesian coordinates satisfies the 
following no n l i near partial differential equation 

y. j. 

(1 - Mo 2 ) ( PxX + <Pyy + Vzz = ^O ^X^XX W 

« 

where U Q and Mq refer to the velocity and Mach number of the undisturbed 
flow, 7 is the ratio of specific heats (7 = 1.4 for air), and x, y, and z 
are Cartesian coordinates. The perturbation velocity vector is given by 
the gradient of cp and has components u, v, and w along the three axes. 

Knowledge of methods for obtaining solutions of equation (l) is 
meager not only because it iB nonlinear, but because it changes type 
(elliptic, parabolic, hyperbolic) . This change of type is an essential 
feature of transonic flow, since subsonic flows are represented by ellip- 
tic equations and supersonic flows by hyperbolic equations. If both 
types of flow occur in a single flow field, it is apparent that the 
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differential equation must change type. In the present case, the type of 
the equation is recognized by the sign of the total coefficient of cp^, 
as follows 

I > 0 elliptic (subsonic) 

(l - Mo 2 )- Mq 2 (7 + 1 ) =~ i =0 parabolic (sonic) ( 2 ) 

< 0 hyperbolic (supersonic) 

• » 

In most of the investigations of two-dimensional transonic flows 
(<p = 0 ), the differential equation is transformed into a linear equa- 

tion of mixed type (Trieomi equation) by the introduction of the hodo- 
graph variables. At the present time, however, no transformation is 
known that achieves a corresponding linearization of the three-dimensional 
equation, and the investigation of other methods of solution thus becomes 
relatively more important. 

Equation (l) is, of course, valid omly in regions where the necessary 
derivatives exist and are continuous. Since these conditions do not hold 
where shock waves occur, and since shock waves are a prominent feature of 
most transonic flows, an additional equation is needed for the transition 
through the shock. The fundamental properties of a shock surface require 
that the normal component of velocity be discontinuous and the tangential 
component, and therefore cp, be continuous. The necessary relation fol- 
lows from the classical expression for the shock polar, which in the small 
disturbance transonic theory is approximated by 

(1 - + (iy a - iVb) + (% ‘ %) 


- Mo' 


, <P X + tPj, 

2 7+ 1 a x b 


(% ■ *0 (3) 


where the subscripts a and b refer to conditions ahead of and behind 
the shock. 

Equations (l) and ( 3 ) are usually developed 
for the case where the coordinate system is , z 

placed so that the x axis is parallel to the J 

undisturbed stream at infinity, but they also > ^ 


a-PPly to the case where the coordinate system 
is rotated slightly. In the present analysis, 
it is convenient to aline the x axis with the 
longitudinal, axis of the wing or body as shown 
in sketch (a) . Such a system is usually 
referred to as the body axes. With these coor- 
dinates, the relation between the total velocity 
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potential ®(x,y,z) and the perturbation velocity potential <p(x,y,z) is 
approximated by 


®(x,y,z) = U Q (x + oz) + cp(x,y,z) 


w 


where a is the angle of attack. 

The expression for the pressure coefficient Cp is not invariant 
with respect to small rotations of the coordinate system. In body axes, 
the proper expression is 

Cp = - ”* (q>x + ctfp z ) - (qy 2 + cp z 2 ) (5) 

u o U Q 

The boundary conditions require that the gradient of the total 
velocity potential evaluated infinitely far from the aircraft be con- 
sistent with the uniform free-stream conditions there and, when evaluated 
normal to and on the surface of the airplane itself, be zero. The condi- 
tion at infinity yields $(«>) = U Q (x + az) or that 


<P(°°) = 0 (6) 

An exception to this statement occurs in the vicinity of the wake at 
great distances behind the wing, but no complication ensues due to the 
relative smallness of this region. The condition at the airplane surface 
results in the relation 


+ni S + n2 ! + “ 3 tb 0 (T) 


where n x , n^, and n 3 are the direction cosines of a normal to the air- 
plane surface with respect to the x, y, and z axes, respectively. This 
relation is too general for the present needs, however, because it applies 
to all shapes, whereas the analysis is to be a small disturbance theory 
that applies only to slender bodies and thin wings. For such configura- 
tions, n x is small nearly everywhere on the surface and will be neglected 
in comparison with either unity or (n^ 2 + n 3 2 ) 1/2 . In this way, equa- 
tion (7) simplifies to 

U 0 (ni + ctn 3 ) + n 2 ^ + n 3 ~ = U 0 (ni + om 3 ) + ^ = 0 (8) 

where n is the normal to the curve bounding a cross section in a plane 
normal to the x axis . 



All of the subsequent analysis proceeds from Green's theorem which 
relates a volume integral over a region V to a surface integral over 
the surface E enclosing V. Green's theorem can be expressed in many 
ways; here it is found convenient to use the forms associated with the 
linear differential equation obtained by equating the left-hand member 
of equation (l) to zero. This results in two different forms of Green's 
theorem, one for Mq < 1, and the other for > 1 and p rom pt s the 
introduction of the following abbreviations " 

0 = (|l - Mo 2 |) l/2 , h = Mo 2 ^i (9) 

u o 

If the undisturbed flow at infinity is subsonic (i.e., Mq < l), 
equation (l) can be rewritten as ~ 


f> a< Pxx + V + * WA* - * 5 (t0 ( 10 > 

and the corresponding expression of Green's theorem is (see, e.g., ref. 21) 

JJJ [*L(Q) - fiL(i|f)]dV = - JJ(y |r - ft |j)dE (ll) 

V S 

where ft and i|r are arbitrary functions and L(fi) is defined as follows 

L(ft) s 0 2 ®xx + %y + ^zz (12) 


and dft/civ is a derivative along the conormal and is defined by 


bSl _ Sft _p Sft dft 

= + (13) 

where n 17 n 2 , n 3 are the direction cosines of the normal to the surface 
drawn into the region V. 

If the undisturbed flow at infinity is supersonic (i.e., Mq > l), 
equation (l) can be written as 
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and the corresponding expression of Green’s theorem is 


JJJ C*L(ft) - flL(ilf)]dV = - JJ(y || - ft H)dE (!5) 

V s 

where the following definitions 


L(0) = *■ P ^xx ^yy ^zz 


and 


apply. 


8ft 

8v 


8ft 

8x 


2 8ft 8ft 

p m + ^ n 2 + ^ n 3 


(16) 


(17) 


DERIVATION OF INTEGRAL EQUATIONS FOR TRANSONIC FLOW 


In this section, integral equations corresponding to the transonic 
differential equation are derived for subsonic and supersonic free-stream 
conditions. One of the principal contributions here evolves from the 
attention given to the shock waves, or discontinuity surfaces, appearing 
in the flow fields. It will appear (see eqs. ( 23 ) and (30)) that the 
perturbation velocity potential can be expressed, for Mq less than or 
greater than 1, as the sum of integrals that show no explicit contribution 
of the shocks. Closer analysis of these integrals reveals, however, that 
discontinuities in velocity can appear and that they automatically satisfy 
the shock-polar relations (see eqs . (3*0 and (37) ) • This section is pref- 
atory to the formulation of the transonic integral equations for the par- 
ticular cases of a slender body of revolution and a thin wing. 


Integral Equation, Mq < 1 


The function in Green’s theorem, equation (ll), is now identified 
with the fundamental solution l/tf of the differential equation L('4 r ) = 0 
and the function ft is replaced by 9 , the perturbation velocity potential 
of the flow field under consideration. From equations (10) and (ll), the 
following relations hold 



- * It £)“- -HJ° l(9)w w 4v 


(18) 
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where 


(T = [(x - x x ) 2 + 0 2 (y - y^ 2 + p 2 (z - z 1 ) 2 ] l/2 (19) 


In these equations the running coordinates in the integrations are x x , 
y x , and z x and cp is to "be calculated at a point P with coordinates 
x, y, and z. 


Equation (l8) is now applied to the inf init e region V surrounding 
the given object to be studied. Same care must be exercised, however, 
in fixing the enclosing surface 2 since Green’s theorem requires that 
singularities and regions of discontinuity must be excluded from the 
domain of integration. It is to- be noted, first, that a vanishes at 
x = x i> y = and z = z x and the effect of the resultant singularity 
can be determined only after the field point is enclosed by a neighboring 
surface and the region V taken external to this surface . Second, since 
shock waves are to be expected wit hin the flow field and discontinuities 
in the perturbation velocity components occur across these waves, the 
boundary of V must also be drawn so as to exclude such discontinuity 
surfaces . 


In sketch (b) , a schematic 
indication of the body and the 
region of ' integration is shown. 

The complete three-dimensional 
extent of the body has not been 
pictured; it suffices, however, to 
state that the surface 2 (shown ' 
dashed) is composed of a sphere of 
large radius which forms the 
external boundary of V, a sphere 
of infinitesimal radius surround- 
ing the field point P, and a final 
surface enveloping the object, its 
wake, and its shock waves. 

If equation (l8) is applied to 

this region and the a priori as sump- 

tion is made that the perturbation 

field attenuates suff iciently fast - Sketch (b) 

with distance to negate the contribution of the surface integral over the 
large sphere in the limit as the radius goes to infini ty, the following 
expression results: 



9(x,y,z) 






dV 


( 20 ) 
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In this equation, the integration region over the surface of the object 
and its wake is denoted by 0 + ¥. The derivatives in the surface inte- 
grals are, in all cases, along lines directed away from the integration 
surface and into the three-dimensional domain V since, as follows from 
equation ( 13 ), the direction numbers have the same sign as the direction 
numbers of the true normal and the Mach number effect is limited to a 
foreshortening of the longitudinal dimension. On the shock surface ~k 
the conormals are directly opposed on the upstream and downstream faces 
Tsq, and Ajj. On the body itself, the conormal derivative can be simplified 
in the manner used in developing the boundary conditions of equations (8); 
that is, from the restrictions imposed on the gradients along the body 
surface, the direction of the conormal becomes effectively that of the 
normal n lying in the crossplanes x x = const. Thus 


dv 




jL_ 

dz 


JL 

dn 


on the surface of the body and wake. 

If the triple (spatial) integral of equation (20) is integrated by 
parts x-wd.se, the resultant form is 


<P(x,y,z) - - 



0+W 


- ! » I; - 


' t ' fe " ni ) - 9 ?] dE - 

^a 

iff &&-!■ ->■■■)" i * 

T'b 


(21) 


Equation (21) is of particular interest because the Integrals over the 
shock surfaces may be shown to vanish. In order to prove this, one notes 
first that the two integrals extend over the same geometric surface but 
that the integrands are evaluated, respectively, on the upstream and down- 
stream faces and, by definition, the directions of the conormals are 
opposed. When the integrands are combined, the total integrand can be 
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expressed as the difference of the two terms, one of which contains the 
factor (<pk - (<p), and the other contains the factor 
A b *a 



( Py 1 n 2 l + ^3 







The first of these factors vanishes by virtue of the fact noted previously 
that the perturbation potential is continuous across a shock surface. The 
second factor can be rewritten in the form 


• %) (nj *b + (%> ■ + (% - w 0 (ns) ^ - 

t(v-v) (ni) *b 

Also, the change in the velocity vector occurring at the shock surface 
must be in a direction normal to the shock. This implies the relations 


(ni) *b 1 ( “ a) Ab = ( “ 3 \ = (% - %) = (% - X) = (% - %) 


Thus, the second factor to be evaluated becomes 



The numerator of this fraction, however, corresponds to the shock-polar 
conditions of equation ( 3 ) and the expression vanishes. 

It finally remains to remark that in the surface integral over the 
body itself , the term nffe-j 2 resulting from the integration by parts 
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is of higher order than the normal derivative of the perturbation poten- 
tial. The term can therefore he neglected and equation (21) becomes 


cp(x,y,z) = - 


JL_ 

kit 




c Rx 1 


0+W 


(s:£> dV 


V 


( 22 ) 


Equation (22) provides another integral expression for the perturbation 
velocity potential, for Mq < 1, in transonic flow theory. The first 
integral on the right is algebraically equivalent to the expression for 
cp(x,y,z) in linearized theory and the spatial integral is a contribution 
brought about by the n onlin ear term of the basic differential equation. 

It is of interest to r emar k that a derivation ignoring the existence of 
the shock waves can also lead to the same form of the equation. In this 
respect the relation is not unlike cases aris ing in linearized supersonic 
theory where it becomes necessary to study the contribution provided by 
the foremost shock wave induced by the body. For the majority of cases 
of practical interest, it can be shown that compensating terms arise and 
that the discontinuity surfaces are taken care of by a formal develop- 
ment that ignores the existence of these surfaces (see, e.g., refs. 22 
and 23) • It is not possible, however, to ignore so completely the exist- 
ence of the discontinuity surface in transonic flow and, as will be seen 
in the later developments, equation (20) is, for certain purposes, pref- 
erable to equation (22) . 


Equations (20) and (22) will now be written, for purposes of refer- 
ence, in the following final forms 


- <H(w) + fcfjjk Qk ?) w 

V 

- %<■***) - 1 s JF i dV 


( 23 a) 


(23b) 


■ - £iT A (M) * 32 - fejQTCsr dV (23o) 

v 

= ^' z) -*iiT A ^)“ dE -Ksiry r¥> dV < 23d > 
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In each of the above equations, has the analytic representation it 
has in linear theory. The first two relations are obvious repetitions 
of equation (22) . The two latter relations are transcriptions of equa- 
tion (20) where the notation A(dcp/dv) = (dcp/dv)^ + (Scp/Sv)^ has been 

introduced, the continuity of (p at the shock surface has been used, and 
where in equation (23d) the variable 


u = sinh -1 


x - x-. 


U .1, LJLl I _ . _ 

pSy-y x ) 2 + (z- zj 2 ] y 

(x-xj |x-x x | + {(x-x x ) 2 + p 2 [(y-y x ) 2 + (z-z x ) 2 j)- 
l*-*xl n Pt(y- y x ) 2 + (z- z x ) 2 ] 1/2 


is employed to express the integral equation in a form that will be of 
value in establishing a reduction to the case of sonic flow. 

The longitudinal perturbation velocity is given by the x-wise 
derivative of any of equations (23) . Consider, for example, equa- 
tion (23a) . After first isolating the singularity at the field point 
by introducing the limits x x = x ± e, one has 

u (x,y,z) 

■ s ^ x ’ y ’ z) + 6 iSo t hJJ iy ^U I 


I \ fe x a (axj. CKI 


= u L (x,y,z)+ g lto o ^JJ 1 |cp Xi 2 (x;y 1 ,z 1 ) 


edy-Ldzx 

[ 6 2 + p2(y-y x ) 2 + p 2 (z- z x ) 2 ] S/S 




Tn the limi t as e— >0 the influence function in the integrand of the 
double integral is effectively a two-dimensional pulse function at the 
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point y x = y, z x = z and of strehgth 2j t/p 2 . The expression for u 
then "becomes 

u (w , . i -hMSOk *> <*> 

V 

A detailed account of -the application of the two-dimensional form of this 
equation to the calculation of airfoil pressure distributions has teen 
given by Spreiter and Alksne (ref . 15) . 


Integral Equation, Mq > 1 


Use is now made of Green’s theorem as expressed in equation (15) and 
ft is set equal to the perturbation velocity potential cp(x,y,z) . The 
direct analogue of the derivation in the previous section would require 

1/2 

that i]r be replaced by [(x - x x ) 2 - p 2 (y - y x ) 2 - p 2 (z - z x ) 2 ] but 

this leads to the immediate introduction of a finite-part technique in 
the Integration. For the initial stages of the analysis, will be 
identified with the fundamental solution to of L(ijr) = 0 used by 
Volterra (ref . 24) . From equations (l4) and (15) , the following rela- 
tions hold: 

5 ) d2 - - - jQ3 r “ £ (^~) df (25) 

2 V V 

where 



_ , x - x, 

to = cosh -1 y— 

eUy-yJ 2 + (z-sj 2 ] 1 ' 2 

The successful application of 
equation (25) to the transonic prob- 
lem hinges on the proper choice of 
the three-dimensional region V 
and its enclosing surface £. Dis- 
continuities in the velocity compo- 
nents are again to be taken into 
consideration at the shock waves. 
Furthermore, the fundamental solu- 
tion of Volterra becomes infinite 
at y x = y, z x = z, that is, every- 
where along' the line passing through - 
the field point P and parallel to 
the x axis. Sketch (c) indicates 


Sketch (c) 
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the disturbance field of the object as well as V and £ (shown with 
dashed lines). The bow shock fixes the foremost extent to the disturb- 
ance field and 2 lies adjacent to it and other possible shock surfaces 
as well as the surface of the object and its wake. The downstream limits 
of the region V are fixed by the forecone with vertex at P and deter- 
mined explicitly by the relation 

(x - X.J = 3[(y - y x ) s + (z - z^.) 2 ]^ 2 (26) 

The inner boundary of V is the cylindrical surface of infinitesimal 
radius given by the relation 


(y - yi ) 2 + (z - zj 2 = e 2 

The cono rmal derivative is defined by equation (17) ; on the infinites imal 
cylinder its direction is parallel to that of the normal to the surface, 
and on the forecone from P the conormal is directed along the surface 
itself. Formal analysis yields the expression 



where integrals over the surface of the body and wake are denoted by r, 
over the two sides of the shock surfaces by Aa and Ab, and over the 
enclosed volume by V. In each case, only that portion of the surface 
or volume lying within the forecone of P is included in the integrals. 
The surface integrals over the forecone itself vanish because w and 
3 ui/dv are zero. It should be noted that the forecone is that of lin- 
earized theory and has no relationship to the region of dependence in 
the actual flow field. 

Integration by parts, in the last integral, leads to the relation 

<P(*,y,z) - - £ | -JJ]* (jg - | - ft S ] a5 ^ ’ 

T 

b h (If + JD [■ d ' I *>) ' * ^ ■]“* 

^a ^b 

b bJFi ^ (si°) d? 

1 


(28) 
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Equation ( 28 ) is the form, for Mq^ 1 , analogous to equation ( 21 ), for 
Mo < l, and on the tody and wake surfaces involves the approximation 


A 

dv dy 


+ 


dz 


_d_ 

dn 


It is not difficult to show, from the shock-wave relation of equa- 
tion (3), that the combined integrals over the surfaces 7 ^ and Tq, 
vanish. The perturbation velocity potential can, therefore, be given 
alternatively as 


<p(x,y,s) - - Jj § ‘ <p lD aZ:+ at teJIfe *** CO 


d V 


T v < 29 ) 

Equations (27) and (29) may now be written in the various forms 


<p(x,y,z) - it(*.y.*) + S’ W* CO® 


(30a) 


" ' &t IM **** 5 dV 


(30b) 




(30c) 


= <P L (x,y,z) - i JJa (g?) i az ±JJTQ~ %“) j dv (30a) 


where cp^ has the same analy tic form as in linearized supersonic theory 
and use has been made of the relations 


oj = cosh -1 


x - x x 


pUy-yJ 2 + (z-zj) z ] x/z 


do) 


-1 


-1 


[(x-x x ) 2 - 0 a (y-yi) 2 - P^z-Zi) 2 ] 


2,X/2 


( 31 ) 
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Comparison of equations (23) and (30) shows once more a difficulty 
that appears in linearized analysis of subsonic and supersonic flow, 
namely, that complete parallelism between the formulas is not achieved 
directly. This parallelism can only be established after interchanging 
the order of integration and differentiation and, because of the singu- 
larities involved, it becomes necessary to introduce the concept of 
finite-part integration. Furthermore, it is well known that the resulting 
multiple integrals can no longer be written in a unique form (see, e.g., 
ref. 25) but must be expressed differently, depending on the order in 
which the integrations are to be performed. Bfo attempt will be made to 
develop these ideas further at the present time . 


Calculation of Conditions on Shock Surfaces 


It is of same interest to study equations (23a), (23b), (30a), 
and (30b) as the field point approaches a discontinuity surface and to 
discover the mechanism by mean a of which these basic equations furnish 
the velocity jumps associated with the shock waves in the field. To this 
end, consider first the case Mq > 1. 

Sketch (d) shows the geometry of tiie 
problem. The bow wave induced by an 
aerodynamic shape is indicated and the 
point P, at which conditions are to 
be calculated, is chosen arbitrarily 
close to the rear surface of the wave. 

The surface of the wave can be replaced 
locally by a planar element and a new 
coordinate system g, tj, £ introduced 
with the origin fixed at the intersec- 
tion of the line y x = y, z x = z and the 
bow wave. Point P then has the coor- 
dinates g, 0, 0 and the planar surface 

is given by the linear relation Sketch (d) 



ag x + bq x + c£ x = 0 


(32) 


Since the bow wave is situated upstream of the linearized disturb- 
ance field, is zero and, from equation (30a), the perturbation 

velocity is 




(33) 
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where to = cosh _1 [(g - iJ/pOlj® + £ X Z ) 1,Z ]- By virtue of the field 
point’s nearness to the how wave , the term u 2 /2 in the integrand is 
assumed a constant and one then gets 


aYo n Xo 

a,,/ 

05 <7. V Y OS l 


where 


Xx = - (htix + eS-J/a, = i - P(ti x s + ^ 1 2 ) 1 ^ 2 

hCai + c^) ± ag[a 2 S 2 + 2ac§£ x - (a 2 g 2 - h 2 - c 2 )^ 2 ] 


2 7 i/2 


*i>*2 = 


a 2 p 2 - h 2 


ai[c ± (a 2 p 2 _ *2)1/2] 

Z 1* Z 2 “ 

a 2 p 2 - h 2 - c 2 


Integration and differentiation yields 


u(£,0,0) 


ku 2 a 2 

2 -a 2 0 2 + h 2 + c 2 


m 


It remains to show that this result, derived from the integral 
expression for the perturbation velocity potential, is consistent with 
the result one would get from the shock-polar relation of equation (3) • 
At the downstream face of the how wave, equation (3) becomes 


-P 2 ^ 2 + v b 2 + V 2 = ku b 3 /2 


m-w* incremental velocity vector occurring at the shock surface is, 
however, normal to the shock surface and this yields the relations 
u b : v b : w b = a:b:c. Substitution into the shock polar relation gives 


fa*b 2 ( u b 2 'N = ^b 2 a 2 

2 W - v* 2 - Wb 2 / 2 _ a 2 p2 + b 2 + c 2 


(35) 


1 
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in agre ement with equation ( 3*0 • It therefore follows that the integral 
expression, for Mq > 1, will adapt itself on the shock surface to any 
how wave consistent with given body geometry. This result can also be 
extended, to include any shock wave in the flow field. 

An analogous procedure follows for the case Mq < 1 . Let the shock 
surface in the vicinity of the point P at i , 0,0 be given by equa- 
tion (32) and assume that u 2 is composed of a continuous part and a 
discontinuous part that has the constant value u a 2 ahead, of and u^ 2 
behind the shock. Equation ( 23 a) then yields 


lim [u(S-,0,0) - u(i+,0,0) ] 
§ — >0 


JT 


. d£ ndTj 




[(a6+bi| x + o5x)‘ + > a 


-X& 


(36) 


where the double integration extends over the region of discontinuity. 

If the differentiation with respect to £ is now carried within the 
integral signs and £ allowed to approach zero, the value of the inte- 
gral becomes independent of the original limits of integration. In this 
way one gets 


Ug. - u b = lim [u(S-,0,0) - u(£+,0,0) ] = 
(• — 


ka 2 


a 2 0 2 + b 2 + c 2 


u a 2 ~ u b 2 
2 


( 37 ) 

It can be shown, as previously, that equation ( 37 ) agrees with the result 
given by equation (3), the shock polar condition, for Mq < 1. 


REDUCTION TO SONIC PLOW THEORT 


Tn this section, the previously determined equations will be studied 
in the limit as sonic, free -stream speed is reached. The integral rela- 
tions then ass ume forms that correspond to the nonlinear differential 
equation when 0 = 0. 

Integrated Strengths of External Sources 

It is proposed here to determine a relation that will prove useful 
in the following section in connection with the reduction of the integral 
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equations to the special forms appropriate for Mq = 1 . This relation 
will "be recognized subsequently as connecting the integrated strengths 
of the exterior corrective sources in the crossplane x = x Q = const. 
and the rate of change of "body cross-section area. As a means to this 
end, equation (10) is written in the form 


k HSfensI = f (38) 

3 x 2 dx 


where V 2 is- the two-dimensional Laplacian operator in the transverse 
plane. Each term is then integrated over the entire x = x Q plane 
external to the body. The double integral involving V 2 ? can be par- 
tially integrated and converted into a line integral by application' of 
Green's theorem for a plane 



( 39 ) 


where the line or curvilinear integral extends around the curves C 
enclosing the region of integration of the double integral. For cases 
in which the plane x = x Q does not intersect any shock waves, the 
region of integration can be taken at once as the entire x = x Q plane 
exterior to the body. If the assumption is made, as in linear subsonic 
theory, that the normal gradient of 9 attenuates with lateral distance 
sufficiently fast to suppress the contribution of the curvilinear inte- 
gral along the outer boundary, the boundary conditions of equation (8) 
permit one to equate the line integral along C to -U 0 s'(x 0 ) where 
S'(x 0 ) denotes the longitudinal gradient of body cross-section area. 

The integrated form of equation (38) thus becomes 


*F(s; 4)„ u M I) ^ az - D ° s ' (x °> (l,0) 
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Tn the more general case, however, in which the plane x = x Q intersects 
a shock wave as illustrated in sketch (e) , discontlmii ties occur which 



Sketch (e) 


require that the integration region must he divided into two parts, one 
lyin g between the body and the shock wave and the other extending beyond 
the shock wave to infinity. Application of Green 1 s theorem to each 
region and addition of the separate contributions results in additi o na l 
1 lnp» integrals carried around the two sides of the shock surface. These 
two line integrals can be combined into a single line integral, in which 
case equation (4o) can be written as follows 


*JF(k 


dy dz + 






bu 

& 


x=x^ 


dy dz - UqS ' (x Q ) 

(M) 


where A(dq>/dn) = dcp a /c)n - dcp^/chi and where in the single line integral 
the integration extends around the curve described by the intersection 
of the shock wave and the x = x Q plane, and the normal n is taken as 
directed away from the body. 

In the subsequent work, attention is to be directed toward results 
at Mq = 1. We assume here that in the limit as 0~>O the first term 
in the right-hand member of equation (4l) will vanish and one then gets 


■JJ(k t) 


dy dz + 


x=x Q 



U 0 S'(x 0 ) 


( 42 ) 


It i-Till become evident in the discussion contained in the following 
section that the left side of equation (42) represents the integrated 
str eng ths of the exterior corrective sources in the cross-plane 
x = x Q = const, in the limit as Mq approaches 1. It follows that under 
conditions corresponding to sonic flight speed, the total source strength 
in any transverse plane is zero: the sum of the sources within the body 
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or wing (sources appearing in the term cp^) is of equal magnitude jjut 
opposite sign to the corrective sources required "by the n online ar term 
in the differential equation. 

Equation (42) allows one to make some conclusions about the lateral 
attenuation of the x-wise gradient of u 2 . Consider the double inte- 
gral as written in terms of polar coordinates 

f ****•>■ 

Since the definite integral must converge, it follows that the integrand 
attenuates faster than l/r^ and if the assumption is made that the 
integrand has a purely algebraic character, one concludes that for large 


5 u 2 _ 1 

^ 2 ri 2+N 


(43) 


where N is some positive constant. 

A check on equation (43) is provided by the work of Guderley and 
Yoshihara (ref. 2 6 ) on axially symmetric flow at sonic speed. In that 
analysis, for large r x . 


_d_ uf; „ 1 

Bx 2 ' 16 /T 
r i 

and this is in agreement with equation (43) when N = 2 / 7 . The same 
reference also gives 

dq> _ 1 

> 9/7 

ar x r x 

which serves to substantiate the assumption made earlier that the curvi- 
linear integral along the outer boundary can be neglected in equation (4o). 


Integral Equation for Slender Bodies, Mq = 1 


In the reduction of the integral equations to the case p = 0, 
methods analogous to those employed in references 25 and 27 will be used. 
Attention will be confined here to field points at a finite distance 
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from the body so that in the limit as 0 approaches zero, the term 0r 
can be assumed to approach zero uniformly. As in conventional slender - 
body theory for linearized flow, the longitudinal distribution of cross- 
section area S(x) is assumed to possess a continuous x-irf.se derivative. 
The method of reduction can be exhibited in a sufficiently general form 
if a lifting body of revolution is considered. When Mq < 1, the pertur- 
bation potential for a body of revolution cpg follows from equations (23) 
and for the purposes at hand the form (23d) is preferable. For suffi- 
ciently slender and smooth bodies the term cp^ can be expressed in terms 
of a rectilinear source and doublet distribution and the corrective source 
distribution appearing in the triple integral, then extends over all space 
external to the x axis. If k(x) is the x-wise distribution of doublet 
strength, equation (23d) becomes 


< P B (x,r,0) 


Do r l s'fa)^ 

“■"4 r a ] 1/2 


sin 0 d f (x-acjXxJdXj, 

4ltr dxJ r „ \2 . ft 2 21 1 / 2 
o Lx -xj + 0 r J 


_1_ 

4it 


I M & 


(x-x x ) 
|x-x, I 


In 


lx— xj + [(x-x.,^) 2 + 0 a p TT 2 3 1 / 2 


dS 


% 

_k_ my a ( x - x x) 

4it dxJJJ \&Xi 2 ) Ix-XjJ 

V B 


In 


1 x-x ± | + [ (x-x 1 ) 2 + 0 2 P I x J l/2 

PPl! 


r x dr i dx x de x 


w 


Tn equation (44), cylindrical coordinates x, r, 0 are used where 
r 2 = y 2 + z 2 , 0 = tan -1 (z/y) and the notation 


p-j-3- = [r 2 + r x 2 - 2rr 1 cos(0-0 1 )] 


1/ 2 


has been introduced. 
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Under the imposed conditions, the first two terms in the right 
member of equation (Mv) reduce, as shown for example in reference 27, to 
approximately 


Ho 

2jt 


S*(x) Znr 


tc(x)sin 0 
2jtr 


Ho 1. 

bn dx 



In 



dx 


1 


It remains to attempt a corresponding modification of the two remaining 
integrals. Consider, next, the triple integral. Since the integration 
extends over all points in space , one encounters a nonuniformity of con- 
vergence in the logarithmic influence function when (3 becomes vanish- 
ingly small and r x becomes infinitely large. From equation (43), 
however, it is known that du 2 /dx attenuates rapidly with increasing r x 
and the resultant error in miscalculating the effect of the logarithm 
for large r x is thereby reduced. The triple integral is therefore 
approximated by 


JL iL /7T T- — 

W dxJJj [dx- , 2 

V 


x-x, 2|x-x x | 

- ■■■ • In dx-Lr-Ldr-LdS-L 

x - x-J pp^ 


This term can be rewritten as 




Consider, finally, the integral over the shock surface. For p 
near zero, one gets the expression 




X - X-l 

|x- Xl | 


In. 


2|x - x x | 


3p. 


dZ 


II 


which becomes 


1 r dq>R / 


£s/t lx - xilaxi / a 
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where dcf^ represents an element of arc on the curve determined by the 
intersection of the shock surface and a plane normal to the x axis. 
The normal n, as in equation (42) , is directed away from the body. 

/ 

The reduced form of equation (44) , for the case Mq = 1, is given 
by the sum of the above three expressions. In the forms given, an 
apparent dependence on 0 remains in the expression 


.11 r . 

4* dxj 


x-x. 2 x-x 

: — 

X-X. 


— dx x U 0 s' (xi) +k JJ 3l ^r 1 dr 1 as 1 + 


T Bci> B 


The bracketed term vanishes by virtue of equation (42), however, and the 
dependence on 0 disappears. The sonic form of the integral equation 
for the ''lifting body of revolution is, therefore. 


< P B (x,r,0) 


u o s ’( x ) ln r _ ic(x)sln 6 + i_ 
2« 2itr 2 it 



Inpj^a-x 


+ 


2it 



In p-g-r^d^de;,. 


(45) 


Starting under the assumption that Mq > 1 and using equation (30d) , 
one can derive the same result, the only essential difference arising 
from the fact that in the limiting process care must be taken to restrict 
the disturbance region V to that portion of space within the Mach fore- 
cone from the field point at x,y,z. 

Equation (45) expresses the sonic equation in the form 

1 P 

<Pg(x,r,e) = cp 2B (x;r,0) + — J A— In p^dcr^ + 

* 

/ 

izJKk ¥) f w lr A (1|6) 

where cp 2T .(x;r,0) is the harmonic potential for the body of revolution 

in transverse planes. As expressed, this sonic form of the integral 
equation for the perturbation potential is identical to the integral 
equation corresponding to the transonic differential equation at 0 = 0; 
that is, it gives, for a flow field in which shock waves may possibly 
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occur, the integral equation corresponding to the partial differential 
equation 

m + m = 7 -- cp q> (47a) 

+yy +zz u Q xxx 

and admits discontinuity surfaces for which the difference relation 



7 + 1 
Uo 2 




is satisfied. The di rect derivation of equation (46) would follow from 
an application of Green’s theorem in the transverse variables to equa- 
tion (47a) without the introduction of the slender -body assumptions. 


The interpretation of equation (42) in terms of net source strengths 
is now apparent. Each of the three terms in the right member of equa- 
tion (46) provides two-dimensional sources in each transverse plane 

= const, and the perturbation flow field is simulated by the combined 
effect of these sources and in lifting cases, a doublet term. (The 
doublet term is of no concern in the present discussion since its net 
source strength is zero.) The first term, q> 2 , contains a source on the 

x axis with strength fixed by the gradient o? area, S ' (x) ; the second 
term represents a curvilinear distribution of sources around the shock 
wave with strengths fixed by A(S<p/^n) ; the last term represents a planar 
distribution of sources with strengths determined by the nonlinear term 
(k/2)(d/tfcc)cp x 2 . Equation (42) thus states that at Mq = 1 the combined 
source strength must vanish in each transverse plane. 

Equation (46) corresponds, for the body of revolution at Mq = 1, 
to equations ( 23 c) and ( 30 d) in that it contains an explicit contribution 
from the shock wave and is expressed in terms of the basic singularities 
of the differential equation. A form analogous to equations (23b) 
and (30b) can also be derived as follows. Let r = R(x,0) be the equa- 
tion of the shock surface. The relation 


p2lt poo g 

a / “T ln 


ps. 1 t 

W 

X 

V. 

as 

p 

~5 

8 

1 

/ flfll 

J + J 

o 

L° J R(x, 01 )J 


— - ■ In p r,ar, 
2 K II 1 1 


p2lt n® 2 

■ s J 401 J t ^ In p n riari “ i - 


k r 2 *./ „ x dR(x,0j f U B 2 \ 

si 


de. 


r 1 =R(x,0 1 ) 
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then holds. Substitution into equation (46) yields 

Inpuddx + 


hSJ ^- to Pn r i dr i as i 

and through use of equation (47b) one gets 



x ' r ’ 8 > = 


m 


The above results have been worked out in some detail for the body 
of revolution. The sonic equations for other shapes follow similarly 
through a reduction of the general transonic integral equations or can 
be expressed directly through consideration of the sonic differential 
equation. The final equations for both cases appear as follows: 

<P(x,y,z) = <P 2 (x;y,z) + 1L.A.JJ ^ In pj^ x dr x dfi x (49a) 

Prft + h Jf^ IT ln enVM®! 

(49b) 


= <P 2 (x;y,z) + — 

" 2 3t <7 


where <p 2 (x;y,z) is a two-dimensional harmonic function which, for the 
body of revolution and for the planar wing is, respectively. 


r2 B 2* 


UqS'(x) K(x)sin 9 

In r - 


2jtr 


(50a) 


>s 2 (x) 


/ ° 2 \ a / 1/2 

Aw w (x; yi )Zn[(y - V X ) Z + z 2 ] cty^ + 


-Si(x) 
ns 2 (x) 

£ J ' , . 

-Si(x) 


z dy x 


[(y-y x ) + z 2 ] 


(50b) 
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In the latter expression, 

Awfoy-J = w(x; yi ) z=o+ - v(x; y;L ) z=o _ 

^p(x; yi ) = 'P(x;yi) z=0+ - 9(x; yi ) z=0 _ 

and the lateral boundaries of the wing plan form are fixed by -s x (x) 
and s 2 (x) . 

SLENDER-WING THEORY IN LINEARIZED FLOW 


In the preceding development, the integral relation for the pertur- 
bation potential in sonic flow has been expressed in a form that follows 
from an application to equation (47a) of Green’s theorem in transverse 
planes. The determination of a solution thus depends to a large extent 
on the evaluation of the effect of the two-dimensional singularities that 
are placed throughout the exterior portion of the flow field. Examples 
of a direct attack on a similar problem are to be found in the calcula- 
tions of two-dimensional transonic flows by Oswatitsch, Gullstrand, and 
Spreiter and Alksne, references 9 through 15- In the present report, an 
indirect attack is to be made, following the ideas of Whitcomb and 
Oswatitsch (refs. 6 and 7) by relating the solution for a slender win g 
to that of a body of revolution. The analysis will show that once one 
establishes the details of the crossflow potential fields associated with 
a wing and its related body of revolution, the residual disturbance fields 
near the two bodies are the same to -a certain order of accuracy in terms 
of the slenderness ratio. The mechanics of such an approach can, in 
fact, be observed in linearized wing theory and such a development will 
be given in this section as a prelude to the subsequent sonic theory. 
Attention will be limited to the subsonic case and, as an added simplifi- 
cation in the analysis, only wings possessing lateral symmetry will be 
considered although such a restriction is not essential. 


Analysis 


The linearized equation for subsonic potential flow is 



+ <P, 


zz 



(51) 
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and if Green’s theorem is applied, formal manipulation leads to the 
following integral relation for the potential cp^ of a planar wing 


9 w (x,y,z) = <P 2w (x;y,z) 



x,r x ,ej)ln [r2 + r x 2 - 


2rr 1 cos(6 - 0 -l)] 1 ^ 2 


ridTidSx 


(52) 


No integrals along possible discontinuity surfaces are necessary since 
shock waves do not appear in linearized subsonic flow theory. Equa- 
tion ( 52 ) is linear in cp and can be separated into additive expressions 
contributing to qfy the potential associated with the thickne ss dis- 
tribution, and cp^ ’ , the potential associated with camber and angle of 
attack. In this way one gets for the perturbation velocity components 

\t 53 u 2w,t ■ £ ^l^+ri 2 -2rr lC os(e -dj] 1 ' 2 x x &r x afi x 

(53a) 

\,t " W % t ~ t ^tr 2+ r l2 -2rr 1 cos(0-e 1 )] l/2 r.dr.de, 

(53b) 

1 

and 

a = ^a, " ^ Jf **i> 2+ r i 2 ~ 2rr lC os (0 - 0 X ) ] l/2 ndridflj. 

(5M . 

w ¥,a = w 2 Wja - fj w , a ^[r 2 + r x 2 - 2rr lC os(0 - 0 X ) f /2 ndr^ 

(5*n>) 


The corrective integrals in equations (53) and (54) obviously do not 
modify the area distribution in the thickness case by virtue of the ver- 
tical symmetry in the flow field nor the load distribution in the camber 
case by virtue of the vertical asymmetry in the field. It follows that 


Aw, 


V,t " Aw =W,t ' **,0. " 


(55) 


where the delta notation denotes the increment in the function in passing 
through the plane z = 0, that is, the difference between the values on 
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the upper and lower surfaces of the wing. As a consequence, cp^ can 
"be expressed in the form 


^(xjyjz) 

Uo 


— r - — — In^r 2 + y x 2 - 2ry.,cos 0) 7 dy + 

2it J U Q 

-s(x) 0 


_1_ r s ( x ) r sin 8 ^ 

2lt -s(x) u o r 2 + y.^ - 2ry.,cos 0 


(56) 


If exact conditions on the wing surface are to he sought in linear- 
ized theory, equations (53a) and (53b) must also he satisfied. For 
example, in the direct case of given thickness, equation ( 5 6 ) predicts 
Ua and equation (53b) is then used to determine the exact streamwise 

W j”C 

velocity component as affected hy the external-source integral. In the 
direct case of given loading, equation ( 56 ) predicts and equa- 

tion (5^h) is then used to calculate the true wing camber, modification 
of w 2 ^ being produced hy the integral term. The difficulties of such 

calculations are so disproportionate to those of solving the linearized 
equation hy standard methods that they appear to add needless complies - 
tions to a relatively simple problem. In transonic theory, however, the 
right-hand member of equation ( 51 ) is replaced hy a nonlinear term and, 
in the absence of more obvious methods of attack, the difficulties 
involved in such an approach become less of a deterrent; in particular, 
the integral forms of the corrective terms are of added interest since 
they are" suited to approximations. The details of such an approach will 
not he considered further at the present time since for slender wings 
the use of a related body of revolution yields information of a sufficient 
order of exactness 'in both linear and nonlinear theory. 

In the following work, a complete knowledge of the solution for a 
body of revolution will be assumed known. Thus, if S(x) is the cross- 
sectional area of the body and the distribution of lift is fixed by k(x) , 
doublet strength per unit of length, the linearized solution for the body 
of revolution is 


( P B (x,r,0) 2 ^ p 1 S^xJdxj, 

U o [(x-x 1 ) 2 +P 2 r 2 ] 1/2 


P 2 r sin 0 p K(x 1 )dx 1 

^ [(X-X^+P 2 !' 2 ] 372 


( 57 ) 
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and the perturbation potential in the transverse plane is 


<p 2B (x;r,e) s’(x)Znr ic(x)sin 6 

( 58 ) 

U c 2it 2itr 

The integral relations for the wing and the body combine to give 

% ~ <P B = % " % " JJ (‘Pw " ^xx ln[ra+ r x 2 - 2 rr lC os(e - 63 .) ] X/ “r x dr x a 0 x 

(59) 

Equation (59) is exact, subject only to the restrictions of first-order 
perturbation theory, but with added restrictions on the geometry and 
loading it is possible to show that the magnitude of the final integral 
is negligible to a certain order of accuracy. We now assume the wing is 
slender, that is, s(x) is small in comparison with over-all wing length. 
Let, furthermore, the wing and body be of equal length and have identical 
longitudinal distributions of cross-sectional area. This implies 


S 



Aw w (x,y 1 )dy 1 

Uo 


( 60 ) 


and establishes the condition that the distribution of two-dimensional 
source strength in equation ( 58 ) is equal to the strengths integrated in 
the transverse plane of the sources appearing in equation ( 56 ) . It will 
also be convenient to equate in the same manner the doublet strengths in 
those two equations and one is led to the relation 


K 


w - -/ 


s ( x ) ^p 2w (x,y 1 )dy 1 


-s(x) 


U r 


(61) 


The first objective will be to show that for field points in the 
vicinity of the slender wing the first two terms in the right member of 
equation (59) differ from the left member by an amount that is of higher 
order in s/z . The evaluation of the error term can be perf ormed by 
an iterative process in which the first step starts with the approximation 


<P, 


W 


9 - 9 +9. 


-w 


■B 


B 


( 62 ) 
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Before integration, the integrand in equation (59) will he written 
as the product of Fourier expansions of its terms. For the logarithmic 
term, one has 


1/2 


m=x 


*1 


ln[r 2 + rx 2 - 2rr x cos(0 - S^,) ] =< 


and, similarly. 


V* A* A cos m (® - 0 x) 

lnI ’L W m ' r *2 r 

(63) 


m=i 


r sin 0 


3^ X (^) 8inm0 ' W ^ r 


r 2 + y x 2 ~ 2ry x cos Q 


-< 


ID^l 


00 

£ slnine ’ r Sl yi l 


(64) 


m=i 


If equations ( 63 ) and ( 6 ^) are used, together with equation ( 56 ), and 
conditions of bilateral symmetry are imposed, one gets 


( p 2 w( x Jy> z ) S’ (x) 


U Q l 


2nl 


y a^mCxJcos 2m0 


£-1 ( r /s) 

m=i 


zrn 


+ Y b sm~iM 8 M2m-l)e 
+ L (r/s) 2 " 1 - 1 


m=x 


where 


a2m(x) = ~- 
4jtml 


y' s(x) Aw g (x,y x ) 


-s(x) 


and 


s < r ( 65 ) 


(66a) 




-s(x) 


W 


(66h) 
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Equation (65) holds true "beyond the circle of radius s(x) enclosing the 
transverse section of the wing; within this circle, the expression is 


jfegterg) . lns + f W*) r/r\ m+1 _ J 

U l 2x1 L (2n + l) 2 LW J 

o n=o ' 


cos 2m 9 


Z_. 2m z_i 

m=i n=o 




2m 


1^] + I s±n(2m - 1)e Z W*)[(f) 

m=x n=o 


2n - 2m + 1 


(if 1 ( 1 + 1 XI 

W \2m+ 2n - 1 2m - 2n - V J ' 


r < s 


(67) 


where the coefficients Agp and Bp n are related to the boundary condi- 
tions through the expansions 


Aw a (x,y x ) 

U 0 



jrtAgnCx) 


n=o 


(68a) 


A<p a (x,y x ) 

V 


I-WW (t ) 2 " 1 

n=o 


(68b) 


Once the coefficients agm. bpm-i, Ap n , and B ?n are related in 
magnitude to the geometry of the wing, the size of the integral term in 
equation (59) can be estimated. Since Aw(x,y)/u 0 is proportional to 
t(x)/l, the wing's thickness ratio, it follows from equations (66) and (68) 
that the following order estimates hold 


&2m( x ) = 0(ts/2 2 ), A^x) = 0(ts/l 2 ) 


(69) 
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Let, furthermore, Acp(x,y) /U 0 "be assumed proportional to a(x)s(x) 
where a(x) is a measure of local angle of attack or camber; equa- 
tions (66) and (68) then yield for the remaining coefficients 


^an-iC*)' = O(os/l) , B an (x) = O(os/l) 


Equation (67) can also he written in the simplified form 


m-x 


00 

^ sin(2m - l)e%(x, , 


r < s 


m=x 


with the order estimates 


(70) 


(71) 


F^x, §) = 0(ts/z 2 ) , Gm^x, §) = 0(ts/l 2 ) , %(x, §) = O(os/z) (72) 

Evaluation of order of error in equation (62) .- Equation (59) is now 
written in the form 


<Pw(x,y,z) = ^(xjy^)- - ({^(xjr,©) + cp^x^,©) + l(x,r,0) ( 73 ) 


where 


.00 _ 2 JT 


I(x,r,«) = -^y J ^(cpjj -cpgjlnt^ + r! 2 - 2rr 1 cos(6 -eiJl^riarid©! 

(74) 


o o 


and an estimation of the order of magnitude of l(x,r,0) is to he mnrip 
for field points in the vicinity of the wing. The approximation of 
equation (62) is to he used together with the given expansions of the 



•fcwn-fl - fmprnH rami- perturbation potentials • Tt guff ices to simplify tlie analysis and. estimate the order 
of the error at r equal to s. Cne then gets 

i(x,s,e) n ii(x,B,e) t i g (x,a,e) 

u 0 i u 0 i v Q i 

pP pS 3 t r „ / „ / V 2111+1 -1 r- 

■ " J COB + Wi( x )fej sin^+UeJ hn r x - 

s ° m»i -* *" 

> j 

i <*>* - £/■/" [® . 4- ?) * 

n=i o o *■ 1 ' ' 


I a) t J amtan-D^x, a)] [la a - J (^)“ ZL &' : , 

m=i el=i 

(75) 

After integration with reapect to 0 X one has 

I( ' * ,6,8) « £ [ f (iP aoa Safi + £®ttl (£)“** aln(2nrt-X)elr,dr. - 

UJ' * O z-> L an Vi/ 2m+i w ■ J - - 

O a Bal w 

"*/ f + rfx, il)lln a r ltol + £ f V R^T ^ + 


I X ~ / J 


CO Z-J L\“/ lim* ~\ ' B / 

0 m=i 


ain(2m - 
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From equations ( 69 ), ( 70 ) , and ( 72 ), the order of magnitude of l(x,s ,9)/U 0 l 
is given by 


l(x,s, 0 ) 
U 0 2 ' 



(76) 


The first step in the iteration fixes the maxim u m value of the error 
inc urred in neglecting the integral term of equation (7*0 • la ‘ftie vicin- 
ity of the wing, therefore, the perturbation potential can be expressed as 
in equation ( 62 ) with an error of the order given in equation ( 76 ) . 

Near the body, a further reduction of the difference - cp 2 ^ is 

possible since the explicit equations for the body of revolution are 
available. This yields the usual form of the slender -wing solution. 

Thus, from equation ( 62 ), 


% ~ <P 2 W + f (x) 


(77) 


where 


f(x) = lim (<Pb “ <P~) — §T f f 

r — >0 ^ /-1 '•v, 


,23 1 

"O O 


(78) 


From equation (57), f(x) becomes, for the subsonic case, 

a 


i _ OO 

f(x) ( X “ X l^ S, ( X 3-) dX 3. J 


In. rx^x^i 


o [(x-xx ) 2 + p 2 rx 2 ] 3/2 


Z 2 

r in uldxx 

83 tSxJ Ix-Xxl V1 L P 2 J 


s . u oi r 

^ 3 t 8 x J 


|x-xxl 

1 s'Cx.na 

Ix-xJ f> 


dx-, 


(79a) 


It is not difficult to show by a similar analysis taking into account 
the possibility of discontinuities in the flow that equation (77) holds 
also for Mq > 1 and that f (x) then has the form 

■■ x 2(x-Xx) 




p 


dx. 


(79b) 
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SLEKDER-raiG THEORY IN SONIC FLOW 


The extension of the foregoing derivation to the nonlinear case 
will he given in this section. The iteration procedure designed to dis- 
cuss the linear problem is effectively the development of an exp ans ion 
in terms of the slenderness parameter of the wing and with appropriate 
restrictions on wing angle of attack can be applied with little 
modification. 


Analysis 


Equation (49b), the relation fundamental to the following discussion, 
provides that the sonic expression for 9 valid in the x = x Q plane 
can be considered to be composed of three terms; cp 2 , a line .integral 
around the possible intersections of the x = x Q plane and a shock wave, 
and a surface integral over the entire portion of the x = x Q plane 
exterior to the body. In many important cases, simplification occurs 
because the line integral introduces no contribution to the values for 9 
in the vicinity of the body. Perhaps the simplest case in which this 
situation develops is that encountered very frequently at Mq = 1 in 
which the shock waves are situated entirely downstream of the most rear- 
ward point. A second case in which the line integral introduces no con- 
tribution occurs when the discontinuity surface is situated in an x = x Q 
plane and is, therefore, essentially a normal shock wave. The disconti- 
nuities associated with the normal shock wave are contained in the con- 
tribution of the double integral. Since most sonic flows about smooth 
wings or bodies probably fall into one of these two cases, attention will 
be confined in the following discussion to those cases in which no con- 
tribution results from the line integral. Thus, if equation (49b) is 
written first for a wing and then again for a body of revolution, and 
the latter is subtracted from the former, the following relation is 
obtained: 

9 w (x,y,z) = 9 (x;y,z) - 9 (x;r,0) + 9g(x,r,0) + J(x,r,0) ( 80 ) 

W -D 

where 

„co _ _ _.l 

J(x',r, 0 ) = —J J | ^ (V 2 - u B 2 Vn U 2 - r i 2 “ 2 rricos (0 - 0 jTrjarja©! 

00 J 

(81) 

The quadratic nature of the integrand in equation (&L), together 
with the additive dependence on thickness and camber in the transverse- 
plane potentials, prompts one to simplify the analysis to cases involving 
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a thin wing of given thickness but limited to an angle of attack or 
camber-length ratio a that is smal 1. in comparison with the thickness - 
length ratio t/l. In this way, sufficient information is retained to 
establish the relationship between the wing and body flow fields for the 
thickness case and at the same time determine a linear dependence on 
angle of attack of the wing loading in the vicinity of a = 0. Under 
these conditions it will be possible to relate the wing flow field to 
that of a body having the same area distribution but not inclined to the 
free stream. Consistent with these conditions, the following equations 
apply: 


92 W " %,t 


+ <P ; 




" 92 B,t 


( 82 ) 


and the perturbation potentials for the wing and body can be expressed 
as 


< Pff * %l,t + ‘Pffjcx. ' 'ffe * <PB,t 


(83) 


where the subscripts t and a identify the contributions attributable 
to thickness and camber. The term u^ 2 - u B 2 in the integrand of equa- 
tion (8l) can now be approximated by 


- V * («aw,t - u =B,t) + - u =B,t) 


(84) 


where the initial assumption 


"v ■ % - v* + (85) 

has been made and higher order dependence on a has been deleted. It 
remains to show, through the evaluation of j(x,r,0), that the assumption 
made in equation ( 85 ) holds. It should be noted that, to the order of 
exactness of this equation, normal, shocks on the wing and body are sit- 
uated at the same longitudinal station. 

The first two terms in the right member of equation (8k) depend 
solely on the thickness distribution of the wing and body, and the two 
remaining terms cont ain the effects attributable to the lift and thickness 
combined. Substitutin g from equation (8k) into equation (8l) , we see 



that the first two terms contribute to J(x,r,0) a function that Is symmetric about z «* 0 and. the two 
latter terms contribute an asymmetric quantity. From equation (80) one then gets 


A %t - = ^ w , a 

These relations are identical to those given in equation (55) and, as a consequence, equation (5 6 ) 
necessarily remains valid along with the expansions given in equations (65) and (67) . 

We now approximate j(x,r,0) at r = s: 
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The orthogonality of the trigonometric terms, together with the fact that Ug ^ has no depend- S" 
ence on 0, permits me to achieve considerable simplification after the integration. The 
resulting expression is 


w A A f 

Ur, l 




(^1/ s) 


m«i 


fri/s) 


sin(2m- 1)0 b^ m . 1 ( x ) 
(2m -1) (r./s)^ 


y---l -**"(■■?)• 
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r 1 dx 1 + 


The convergence of all the terms in the first integral is assured when ug ^ varies as l/r** for 
large r, N being any positive constant; for small r we assume Ug ^ varies directly with cross- ' 
sectional area of the body. The order of magnitude of J(x,b, 0)/U o I 'can then be seen to be given by 


j(x,s,0) 


U Q l 


= ° (^ ta s) 


(87) 


Since it follows from equations (85) and (56) that the magnitude of Vyfu 0 l in the vicinity of the 
wing is 0[(ts/Z 2 )Zn s], the relative error incurred by neglecting j/u 0 Z is 0(ts 3 /z 4 ) . 
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Following the method used in deriving equation (77) , we can achieve 
a final simplification in the wing's perturbation potential. Thus, 

<Ptf * <P 2W + g(x) (88) 

where 

/ \ r°° S u, R 2 ( x ^ r i) 

g(x) = lim - 9 2Bft J = k J ^ “^-g ^ r^cb^ (89) 

r >0 * o 


In the absence of analytical solutions for the body of revolution problem, 
the evaluation of g(x) must be carried out by less direct methods. This 
will be discussed further in the section on applications. 


APPLICATIONS TO SEVERAL PROBLEMS INVOLVING SONIC FLOW 


In the following section, the exploitation of the results for sonic- 
flow conditions will be carried out in some detail. In view of the 
difficulty associated with transonic analysis, it appears likely that the 
equivalence relation of equation (85) will play an important part in the 
interpretation and use of experimental data as well as in purely theoreti- 
cal predictions. The discussion will be concerned principally with appli- 
cations to slender wings and bodies and to the relationships between the 
aerodynamic characteristics of the two configurations. It is obvious 
that the known basic information can be supplied either by theory or by 
experiment and many of the results to be given are written with the idea 
that they can be used in this dual sense. 


Resume of Principal Results of Slender-Body Theory 


It appears worthwhile, before proceeding to the examples, to 
re-examine the problem of transonic flow about slender bodies of arbi- 
trary cross section from a heuristic, although less rigorous, point of 
view. This second approach may be regarded, if one prefers, as a physi- 
cal interpretation of the result given in equation (80) . 
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Consider the case of compressible flow about the slender body ■of 
arbitrary cross section shown in sketch (f ) . The general procedure is 

SLENDER BODY THEORY 



Sketch (f) 


to consider first the complete three -dimens ional problem, and then to 
introduce simplifying assumptions consistent with the restriction of 
slender plan forms. As in the more detailed analysis, it is again assumed 
that any shock waves are situated either entirely upstream of the most 
forward point of the body or entirely downstream of- the most rearward 
point, or are normal shock waves if situated along the length of the body. 
The resulting solution for the perturbation potential can be expressed 
as the sum of four parts. As in the preceding analysis cp^ , (p^ 

and 9 2 -d + 8X6 solutions of the two-dimensional Laplace equation as 


•*B,t 
indicated. 


Thus <p. 


2 W,a 


corresponds to the ■ two-dimensional incompressible- 


flow solution for translation of the cross section, and cp ; 


2 W,t 


to that 


for the growth of the cross section. In addition to satisfying the pre- 
scribed boundary conditions at the body surface, these two terms satisfy 
the requirement that the lateral velocity components (dcp/dy, dcp/dz) vanish 
at infinity. These terms alone do not furnish a satisfactory approxima- 
tion, however, for cases in which S(x) is different from zero because 
<p~ acts like S’(x)Zn r at large r, and hence dcp/cbc is infinite 

w 

at a large lateral distance. This error can be removed, however, by 
subtracting the term <p 2 corresponding to the two-dimensional 

incompressible -flow solution for the growth of a body of revolution 
having the same S(x) as the original body (thereby canceling <p 2 __ . at 

large r) , and adding the three -dimens ional solution cp^ ^ 


W,t 
for flow 
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about the same body of revolution. If <Pg ^ is determined from linear 
theory, the results correspond to the familiar formulas of subsonic and 
supersonic slender-body theory. (This leads, in both cases, to the 
result given in equation ( 62 ) and reduces to equation ( 77 ) where, for 
subsonic flow, f (x) has the form given in equation ( 79 a) and, for super- 
sonic flow, the form given in equation (79b) •) In keeping with the 
previous analysis, this function of x will be denoted by f (x) if it 
is determined from linear theory and by g(x) if determined from tran- 
sonic theory. 

Although the linear-theory approximation is unsatisfactory at Mq= 1, 
the same intuitive procedure can be extended to sonic flow. The desired 
expression follows if 9 b t is determined from the transonic differential 
equation . Thus, as given in equation ( 85 ), one has to a known order of 
accuracy, 

% - + ‘few.t - 'fcB.t + (9°) 

and this result reduces for points near the body to 


" ^aw + g < x) 


(91) 


where 


% - + ’ eW ’ “f 0 (%t - m 


as indicated in sketch (f ) . It is apparent that equation (90) has a dual 
basis for validity and represents either the relation afforded by tran- 
sonic theory for Mq = 1, or that given by linear theory for other Mach 
numbers . The customary restriction to slender wings and bodies must be 
observed in both applications. 

The power and weakness of the present intuitive reasoning is well 
illustrated by the fact that the relation given by equation ( 90 ) is 
found without recourse to the detailed investigation of the earlier 
sections whereas the restriction to small angles of attack that enters 
in the simplification is overlooked. This deficiency stems from the 
fact that it is insufficient to assure that merely the inf ini te veloci- 
ties be removed. Since the space involved is infinite, it is also neces- 
sary that certain integrals of velocity (see eq.. (43)) be finite, and it 
is in connection with the attenuation of the velocities arising from the 
term cp 2 that the deficiency occurs. One could, perhaps, have 
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continued the heuristic reasoning but, once the principal idea has been 
established, formal analysis can be used to establish the restrictions 
and to evaluate the error terms involved. 


Some insight into the validity of the foregoing equations can be 
obtained by examining the numerical solution given by Yoshihara in ref- 
erence 17 for sonic flow about a circular cone-cylinder at zero angle of 
attack. Inasmuch as it was not assumed either explicitly or implicitly 
that the perturbation potential in the vicinity of the body has the form 
indicated by equation (91) (the boundary conditions were satisfied at 
the actual body surface rather than along the body axis) , these results 
are particularly suited in this respect for the investigation of the 
region for which the simplified relation applies. On the other hand, 
the example is not ideal because the sharp comer at the shoulder violates 
the smoothness condition; it is, however, the only case for which a 
theoretical solution is available. Accordingly, sketch (g) has been 



Sketch (g) 

prepared so as to show the variation of cp^ _j./u 0 with r for several 

stations along the length of the cone for the case in which the semiapex 
angle 0 is l/lO radian. Attention is called to the fact that the 
values for <Pg ^ given in this sketch are for a cone of unit length 

whereas the original values given in figure 5 of reference 17 are for a 
cone of length 10. A dotted line is also shown for each station repre- 
senting the values obtained after subtracting cp 2 /U 0 computed by 


?2 B,t = dS lDr = xtor 
U 0 2k dx 100 


( 93 ) 
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from cpig t /u Q at the same point. In order to illustrate further the 
nature of these results, sketch (h) has "been prepared to show the varia- 



Sketch (h) 


tion of - cp2j) ^ with x for various r. The resulting values 

should, according to equation (92), be a function of x n.l nnp for 
s mall r. It can be seen from an examination of these sketches that the 
difference cpg^t - t is indeed very nearly a function of x in 

most of the region for which results are available. Slight deviations 
occur in the immediate vicinity of the nose and at the largest distances 
from the body. The latter departures are so small, however, that it is 
necessary to possess additional information for greater' distance fr om 
the body before one can determine the extent of the region for which 
the g(x) function is applicable. 

Determination of g’(x) in Terms of Pressure Distribution 
on a Konlif ting Body of Revolution 

Although qXg and hence f (x), can be calculated directly by means 

of linear theory for either distinctly subsonic or supersonic flow, gen- 
eral methods are not yet available for the theoretical determination of 
cpg^ in transonic theory. It is evident from its definition, however, 

that g(x) depends only on the longitudinal distribution of cross-section 
area s(x), and that its derivative can be dete rmin ed from simple aero- 
dynamic measurements of the flow about a slender nonlifting body of revo- 
lution having the same S(x) as the given body. From the point of view 
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of applications, nothing is lost in not knowing the actual level of 
g(x), however, since knowledge of its gradient, g’(x), is sufficient 
for the deter mina tion of flow quantities such as velocity and pressure - 
Since the easiest flow quantity to measure is generally the pressure 
distribution on the surface of the body, perhaps the simplest way to 
determine g'(x) is through a relation expressing this quantity in terms 
of the pressure distribution. The necessary relations for the perturba- 
tion velocity potential and the pressure coefficient Cp are 

provided by equations (48) and ( 5 ), which reduce, in the vicinity of a 
slender nonlifting body of revolution, to 


m u o as . /v 

%* = S5 Inr+g(x) 


(5*) 



-2_i.cn + 

U n dx \dx/ 


. 2 . 1 . 

Uq Sx 


*B,t 


S'%) 

4«S 


(95) 


where R(x) represents the radius of the body of revolution and the prime 
denotes differentiation with respect to x. These relations can be com- 
bined to solve for g’(x) in terms of the surface pressures and the 
cross-section area with the following result: 



The cone-cylinder solution of Yoshihara (ref. 17 ) again affords a 
means of illustrating the application of this result at 1^ = 1. Thus, 
sketch (i) has been prepared to illustrate the variations with x of 
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the pressure coefficient 1 on the surface of a cone-cylinder having a 
semiapex angle of l/lO radian, and of g’(x)/U 0 computed therefrom, 
using equation (96) . As in the case of sketches (g) and (h) , the values 
of g*(x) have been converted from those given ori ginally for a cone of 
length 10 to those for a cone of unit length. 

■ It is likewise evident that the function g’(x) can also he deter- 
mined from pressure-distribution data for thin wings in an analogous 
manner, although naturally more geometric quantities are involved in the 
calculation. 


Relation Between Pressure Distributions - 
on Related Wings and Bodies 


Wings and bodies having same longitudinal distribution of cross- 
section area .- Equation ( 90 ) displays the relationship between the per- 
turbation velocity potential cp^ for a thin low-aspect-ratio wing and 
the corresponding potential cpg ^ for a slender nonlifting body of 
revolution having the same longitudinal distribution of cross-section 
area. In most practical applications, however, one is not so much 
interested in relations involving the velocity potential as those involv- 


ing the pressure distributions. The 
cerned with the derivation of such a 
objects Illustrated in sketch (j). 
Both have the same s(x), but the 
first is a nonlifting body of revo- 
lution and the second is a thin 
lifting wing. The relations for the 
potential and pressure coefficient 
for the body of revolution are those 
given in equations (9*0 and (95) of 
the preceding section. The corre- 
sponding relations for qty and Cp__ 
in the vicinity of the wing are w 




following discussion will be con- 
relation. Thus, consider the two 



o _ 2 ^w 

Cp w u Q ax 


(98) 


1 The curve for Cp^ shown in sketch (i) differs from that given 

origi nally in reference if due to the correction of a sign error in the 
quadratic term of the expression for pressure coefficient. For the cone 
having a semiapex angle of l/lO radian, this change has the effect of 
d imini shing the values given originally for Cp™ t on the cone surface 
by a constant amount, namely, 0.02. 
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Since g(x) is the same for "both objects, the desired relation between 
the pressure distributions on the wing and body of revolution can be 
determined by combining equations (95) through (98); thus 



2 2 W S"(x) S S* (x) 

— — _____ + — — in — + — 

U 0 dx 2 it it 4 jtS 


( 99 ) 


It is interesting to note that this relation holds not only for nonlinear 
theory of sonic flow, but also for linearized slender-body theory for 
subsonic and supersonic flow. This follows directly from the fact that 
equation (91) and the associated statement are equally correct in lin ear- 
ized slender-body theory if g(x) is replaced by f(x). 

The term involving 9^. can be considered known inasmuch as it can 
be determined directly using equation (50) or any of several other 
methods (e.g., conformal mapping, etc.) available from classical two- 
dimensional potential theory, or indir ectly if either the linear theory 
or the slender -wing-theory solutions are known for the wing. To illus- 
trate, let the subscript S denote the values indicated by the slender- 
wing-theory solution. Then, for example, if (qfy) is known, 9^ is 
given by ® 


% = (<Ptf) s - fW 


(100) 


where f(x) is provided by equations ( 79 )* Correspondingly, one has 



+ 


— f’(x) 

u 0 


( 101 ) 


If , on the other hand, the linear -theory solution is available, the 
relation 

lim 9 t = lim 9 g (102) 


applies, whence 


_ 89 

— § — - 3 l = lim 
U Q dx Mq->x 



( 103 ) 


where the subscript L refers to values given by linear theory. 
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Equation (99) enables one to caJ.cu3.ate the pressures in the vicinity 
of any thin law-aspect-ratio wing, provided the pressure distribution is 
known on the surface of a' nonlifting body of revolution having the same 
longitudinal distribution of cross-section area S(x). The corresponding 
rule relat ing the pressures on two wings having different cross-section 
shapes but the same S(x) can be easily derived by applying equation (99) 
twice and subtracting so as to eliminate all terms pertinent to the body 
of revolution. 


Wings and bodies having similar longitudinal distribution of cross- 
section area .- It is a simple matter to extend the previous results so 
as to include more general relations which apply to wings and bodies 
having longitudinal distributions of cross-section area that are merely 
proportional. The information needed to achieve this generalization is 
supplied by the transonic similarity rule for slender bodies of revolu- 
tion (ref . 5) • The rule states that at Mq = 1 the pressure distributions 
on two slender bodies of revolution having area distributions given by 



(3-0*0 


are related accor ding to 



3m, I P B,lW * a (x/lf 



l/2 


(105) 


where Sm refers to the maximum cross-section area and the subscripts I 
and II refer to properties associated with the two bodies. If both bodies 
are in air, 7 -j. - 7-j-j, and equation (105) reduces to 



^.u /xN . 1 d2s n(xA) 

Sm,I ^B,lW * d(x/l) 2 



®m,II 


(106) 


If it is desired to determine the pressures Cp^ for a wing having 
an area distribution given by Sjj(x/l) and the pressure distribution is 
known for a body of revolution having an area distribution Sj(x/z) 
proportional to Sjj(x/l), one merely computes Gp^ ^ for a body of 

revolution of area Sjj(x/l) using equation (106) and substitutes the 
result for Cp-g^ in equation (99) • 
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Application to the Calculation of Pressures and Forces 
on Thin Elliptic Cone-Cylinders at Mq = 1 


Thg relations developed in the preceding section will now he applied 
to the calculation of the pressures and forces at Mq = 1 on the conical 
portion of the thin cone-cylinder of elliptic cross-section illustrated 
in sketch (k) . The or dina tes of the upper surface of the cone are given 

*>y 



+ 1/2 

Zu(x,y) = — (“2x2 - y 2 ) (107) 

where m is the tangent of the semiapex angle, 
l is the length of the cone, and t is the 
maximum thickness of the cone. Xt follows that 
the elliptic section in the plane x = x x has 
major and minor semiaxes equal to mx! and tx x /2Z 
respectively. The cross-section area and surface 
slope are, respectively. 


S(Xi) 


irfcmx-i 2 


2Z 


dgu mtx x 

dxi 2Z(m2x x 2- y 2) ;L/2 


(108) 


Pressure distribution on nonlifting cane-cylinders .- From equa- 
tion (50), cp 2w for the symmetrical nonlifting case can be written as 


% = w ° ihf ^ yi ^ n + z2]1/ ^1 ( lQ 9) 

-mx 

which, when evaluated on the wing surface (i.e., z = 0, -mx < y < mx), 
yields 



U 0 tmx mx 

-■ ■■ ■ in — 

21 2 


( 110 ) 


After inserting this relation into equation (99) 
indicated operations, one obtains 


and carrying out the 


°% = Cp B 


mb /. 

" 21 V 


1 + Zn 



(in) 


where Ct> refers to the pressure. distribution on a circu lar cone- 
cylinder having such a semiapex angle 6 that it has the same longitudinal 
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distribution of cross-section area as the elliptic cone-cylinder j thus 



The pressures on such a body can be dete rmine d from those shown graph!- • 
d aily in sketch (i) for Q = 0.10 by application of equation ( 106 ) , 
which reduces to the form 



= 100 6 



- 2G 2 In 100 0 2 


( 112 ) 



As mentioned previously, following equation (99 ) y the difference 
- Cp is the same as given by linearized slender-body theory for 


subsonic or supersonic flow. As a corroboration of this statement, con- 
sider the expression given by slender-body theory for -the supersonic 
pressure on the thin elliptic cone (ref. 25 , p. 257 ) 


~ = _ nrfc 

P SW l 


1 + In 


m(Mo 2 - 1) 


1/2' 


(113) 


and the corresponding expression for the supersonic pressure distribution 
on the slender circular cone (ref. 25, p. 24l) 


Cp = - 

dB 


1 + 2 In 


e(Mp a - 1) 


1/2' 


mt 

21 


1+ 2 In 


rcfc(Mc 2 - 1) 
8l 


-1 1/2 


(n4) 


The difference between equations (113) and (ll4) obviously reduces to 
the form given in equation (ill) . 


The application of the foregoing theory to a specific case will now 
be illustrated by determining the pressure distribution Cp^ at Mo = 1 
on an e llip tic cone having m = l/2 and t/l = 0.06. The first step is 
to calculate the pressure distribution at Mo = 1 on the surface of a 
circular cone-cylinder having a semiapex angle given by 

N e = ' = ( 0 . 015) 1/2 = 0.1225 (115) 

The pressure distribution on the selected elliptic cone-cylinder can then 
be calculated through use of equation (ill) and is 

% - Oj-B - 


(116) 
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The results for both the circular and the elliptic cone-cylinders are 
shown graphic ally in sketch (l) . Note that the pressure distribution 
is independent of y in this case and that a single curve of Cp^ 
versus x/l suffices to define the pressure on the wing. 

PRESSURE DISTRIBUTION ON AN ELLIPTIC CONE-CYLINDER AT Mo=l 




Drag of nonlifting cone-cylinders. - The drag at Mq = 1 of thin 
P ll-l pbv. cnnR-cyHnders can be obtained by direct integration of the 
product of the pressure and surface slope and is expressible in the form 


% - B e ♦ 

w 

where the integration is extended over the plan form and D e represents 
t he contribution to the drag that results from a f i n i te leading-edge 
radius of curvature. Since only that portion of qfy denoted by cp 2w 

contributes to D e , this q uant ity can be calculated in the same mann er 
as described in references 25 and 28 for linear theory. Thus, the con- 
tribution per unit of span is , in slender theory , 


dDg 

dy 



jt 


2 


(H8) 
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where r n is the radius of curvature normal to the wing lea ding edge 
and s is the local semi. span. If the ordinate of the wing, in the 
vicinity of the leading edge, is 


= h(s,y) (s - y) 


(H9) 


equation (ll8) becomes 


m * PoV .2 

dy 2 2 


/ \ 2 
, J ds\ 


( 120 ) 


For the thin elliptic cone, z u is given by equation (107) and s(x) = mx, 
hence 

h 2 (s,s) = -S2- 
Pm 2 7. 5 


B e = 2 f <Jx. f rft* 

J dy 4 2 


( 121 ) 


The second term on the right of equation (117) becomes, upon substitution 
of equation (ill) for Cp 

If I 1 [/% f - . f (1 + s)b W ] 

w 0 

- ”3 - f ^ “K 1 + i) 

The drag of the- elliptic cone -cylinder is thus 


% = d b " f £ ^ 2 “ m2t2ln H 


( 122 ) 


Note that the circular cone -cylinder and elliptic cone -cylinder have 
different values of drag at Mo = 1, even though they have the same area 
distribution. As an illustration of the order of magnitude of the quan- 
tities involved in equation (122 ) , the drag at Mq = 1 of a circular 
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cone -cylinder having semi apex angle 0 = 0 . 1225 , as determined by inte- 
gration of the pressure distribution shown in sketch (Z) , is 


(Db) 


9=0 . 1225 


0.00484 


p 

PgUp ^2 
2 


(123) 


■whereas that of an elliptic cone-cylinder having m = l/2 and t/Z = 0.06 
is 

(%) = 0.00383 ^2- Z 2 (124) 

H m = 1/2 2 

t/Z=O.OS 


Thus, although both bodies have the same area distribution, the drag of 
the ellip tic cone-cylinder is less than 80 percent of that of the circular 
cone-cylinder . 

More general results for circular and elliptic cone -cylinders can be 
obtained by combining equations ( 123 ) and (124) -with the transonic simi- 
larity mile for the drag of slender- bodies of revolution. The latter can 
be derived by integration of the corresponding relation for the pressure 
given in equation (106) and -was first given by Oswatitsch and Bemdt in 
reference 5* It states that the drags at Mo = 1 of two bodies of revolu- 
tion having area distributions given by equation (104) are related, 
assuming both bodies are in air, so that 7 j = y TT , according to 




Pq U q 2 S I , 2 (7)-S» 2 (0) 

2 2jt 


Zn 


Sm,I 

s m,U 


(125) 


For slender circular cone -cylinders S= it0 2 x 2 and equation (125) 

reduces to 

"n = (^) [^ + ^ ei* 1 * l “ (?£) ] C 126 * 

which becomes, upon substitution of the values given by equation ( 123 ) 
for Dj and 0 j 

Dg = - *0^1.55 + 2 Zn 0 2 J f ( 127 ) 


The general expressions for and Bg can be compared more readily if 
the relation 0 = (mt/2Z ) 1/2 is introduced to express 0 in terms 
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of m, t, and 1 of the thin elliptic cone having the same area distri- 
bution, thus 

% = " [{ m a t 2 (l.55 + 2 In g) ] (12 Q ) 

Combination of equations (128) and (122) yields the corresponding result 
for the drag at Mo = 1 of thin elliptic cone -cylinders 

% - - *§£ [f “ 2 ^(l-55 + l» S£) ] (129) 

Before leaving the subject of similarity rules, it is of interest 
to note that equation (129) for the pressure drag of thin elliptic cone- 
cylinders is in accord with the transonic similarity rule for the pres- 
sure drag of thin wings of finite span (see, e.g., ref. 4). The latter 
is us ually given in dimensionless form and provides that the pressure 
drag coefficient C D at Mq = 1 of a family of thin nonlifting wings of 
affinely related geometry, plan-form area Sp, thickness ratio r, and 
aspect ratio A satisfy the equation 


C D 


T 


5/3 


= f (Ar l/S ) 


(130) 


where 


CD = 


p o U o 2 ~ 

2 ®P 


(131) 


f indicates a functional dependence. If indeterminate forms that 
arise from the infinite plan-form area of a cone-cylinder are avoided by 
letting Sp represent the plan-form area of only the conical part of 
the body, Sp, A, and r are related to m, t, and l according to 


Sp 


ml 2 , A = 4m, t = j 


(132) 


and equation (129) can be rewritten as follows 



% 

p o u o „ 5/3 

T~v 


-2L At 1/s 
16 


C 1 - 


55 + 3 In 



(133) 


It is evident from this form of the result that Cp/T 5 ^ 3 is a function 
of At 1/3 alone, as required by the similarity rule. 
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Pressures and forces on lifting cone-cylinders . - The relations sum- 
marized in sketch (f ) also permit the calculation of the pressure distri- 
bution on a thin elliptic cone -cylinder when inclined at a small angle of 
attack. To calculate this quantity, one must first have the expression 
for cp^ The necessary expression is well known, however, since the 

problem is equivalent mathematically, for the planar boundary conditions, 
to the boundary- value problem associated with translation of a flat plate 
in a two-dimensional inco mpr essible fluid. Thus qu and cp^ on the 
surface of the thin inclined elliptic cone are W,<x 


a = * U 0 a(m2x2-y2) 


1/2 


( 13 ^) 


U n tmx 1/2 

^ ■ ~ 2 ~ ln — ± U Q a(m 2 x 2 - y 2 ) 


(135) 


where the upper (plus) sign is to be used on the upper surface and the 
lower (minus) sign on the lower surface. After inserting equation (135) 
into equation (99) > one obta ins 

where the convention concerning upper and lower signs still holds. The 
aerodynamic loading, or the difference in pressure between the two sides 
of the wing, is thus 


Po U o 2 / 2 



hrm 2 X 

(m 2 x 2 - y 2 ) 1 ^ 2 


(137) 



A sketch of the load distribution is shown 
in sketch (m) . Integration of the loading 
over the plan form leads to the following 
expression for lift 


% = 


2 

(ancon 2 ! 2 ) 


(138) 


Sketch (m) 
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Although the pressure distribution at Mo = 1 is not the same as 
given by linear theory, it will be recognized that the load distribution 
and lift arise solely from cp 2 ^ and are therefore the same as given 

by linear slender -wing theory. One recognizes , consequently, that the 
lift of any low-aspect-ratio wing having such a plan form that no part 
of the trailing edge extends forward of the station of maximum span is 
given by 


% = ^“(2*a£ 0 2 ) (139) ' 

and the drag due to lift by 

% - % " I L ( l4 °) 

" "a=0 2 

The fraction l/2 enters as a result of suction forces on the leading edge. 
Note that the above statements also imply that all reciprocal and reverse 
flow relations of linear theory are applicable to lifting forces at 1 

on slender wings at small angles of attack. 

As in the case of drag discussed previously, equation (138) for the 
lift of thin low-aspect -ratio wings is compatible with the transonic 
similarity rule which states that the lift-curve slopes of a family of 
thin wings of finite span and affinely related geometry are related 
according to (see, e.g., ref. 4) 


T l/S CL a = f (At i/3 ) (l4l) 

where 


= lim * (142) 

a <*->o ( Po U 0 2 /2)S p a 


Substitution of the geometric relations of equation (132) into equa- 
tion (138) for the lift, yields simply 

T i/3n =£at 1 ' 3 (143) 

"a 2 

which is obviously in accord with the similarity rule. The drag due to 
lift given by equation (l4o) is in corresponding agreement with the . 
appropriate transonic similarity rule. 


T 
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Momentum Analysis of Drag 1 of Slender Bodies at Mq = 1 


The previous example of the thin elliptic cone -cylinder has disclosed 
significant differences in the drag at Mq = 1 of elliptic and circular 
cone -cylinders having the same longitudinal distribution of cross-section 
area. Since this finding is contrary to the often quoted transonic area 
rule. It is of interest to study the sonic drag of a more general class 
of bodies. This will now be done using momentum methods. 


Derivation of general relation for drag .- Consider a surface E 
which encloses a volume containing an aerodynamic body. The vectorial 

force F on the body can be determined by considering the pressures and 
flux of momentum at E . In general there results 


F 


-IT*- 


no -> 


If pv 

E 

(^U Q + vj.dE 


(144) 


where vector notation is used, p and p are the local static pressure 
and density, and V is the local perturbation velocity vector. 

For present purposes, the surface E will 
be taken as shown in sketch (n) . Two parts 
of E denoted as I and II, are plane surfaces 
normal to the x axis and situated upstream 
and downstream of the body. The remaining part 
of E, denoted by 111, is a small circular 
cylinder of radius large enough so that 

the body is entirely enclosed within the cylin- 
drical surface. Since it will be assumed that 
the body is slender and smooth enough that the 
necessary restrictions on cp are satisfied .at 
all stations forward of the base, but that discontinuities in geometry 
or velocity may occur there, the plane surfaces I and II will be placed 
at infinity upstream and at the base of the body, respectively. 



It is sufficient, to the order of accuracy of transonic theory, to 
approximate p and p at points near the body by 


Po °o 

and 

P - Po = - Po |u 0 u + | ^1- Mq 2 )u 2 + v 2 + wsj j- 
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Furthermore, if attention is restricted to the streamvri.se component of 
force, total drag of the enclosed body is given by 


D = — JJ j^Mo 2 “ + v 2 + w 2 J dy dz - p Q JJ uv r R£d0 dx (l45) 

n m 


where v r is the radial component of velocity. This expression holds, 
of course, at Mq = 1 and becomes 


D = JJ (v 2 + w 2 )dydz - p Q JJ uv r %d0 dx (l46) 


II 


in 


An alternative form for equation (l46) which will prove useful can be 
obtained by re pla cing the surface integral over H by a line integral. 
Thus, Green's theorem provides 



(147) 


where C is a curve, situated in plane n, which goes around the wing 
a.nfl also around the control surface Z, du c is an element of C, and n 
is the unit normal drawn into the interior of C, as illustrated in 
sketch (o) . But near the body the relation cp = tp 2 + g(x) holds, and 
therefore within C, hence the equation 
^9 =0 is satisfied in II , 


JJ (v 2 + w2)dy dz = - J 9 ^ dh c 

II c 

and equation (l46) for drag becomes 

B = ~ ~qJ ^ ^ d °c - Po JJ uyrRsie dx 

c m 


(148) 


(149) 




6o 


MCA m 3717 


Relation between drag of wings and bodies having the s»m<^ 
distribution .- Consider sonic flow about two aerodynamic shapes, one a 
thin wing and the other a slender nonlifting body of revolution, having 
the same longitudinal distribution of cross-section area. If equa- 
tions (90) and (91) are substituted into equation (l 46 ) , and if the 
Fourier expansion for cp^. - cp 2 ^ obtained from equations (58) and (65) 

is introduced into the integrals over TTT and the portion of C con- 
tiguous with IH, one has 



JJ n 2 


r r Rj<L0 dx 
B 


III 


( 150 ) 


The contour C is here divided into two parts . The inn er portion that 
immediately surrounds the wing is denoted by C^, whereas the outer 
portion is denoted by C TTT . It follows similarly that a corresponding 
expression can be written for the drag of the body of revolution. Thun 




db. 



JN 


T r R v dj 0 dx 
r B S 


III 


( 151 ) 


where Cg refers to an integration contour drawn around the cross section 
of the body. If the exterior portions of the control surface E are 
now selected the same for the two cases (i.e., surface III is the same 
for the wing and body), it follows immediately by subtraction that 


% - % 




(152) 


The integrals over C™ can be divided into two parts after rec alling that 


'Ml 

the analysis applies for sma.1 1 angles of attack and that 


cp> can be 
W 



NA.CA TUT 3717 


61 


written as the sum of cp 2 and cp 2 , provided the thickness does not 

W,t W,oc. 


vanish. Substitution of this relation into equation ( 152 ) yields 



where the integrals involvin g the cross-product terms are absent since 
cp 2 ^ ^ and its normal derivative along Cfy- are even functions of z, 

whereas a and its normal derivative are odd functions of z. 

Since cp 2 ^. is proportional, to a, it is evident that the first 

integral of equation (153) provides a contribution to the drag which is 
proportional to the square of the angle of attack. This quantity is 
exactly the vortex drag and is represented by the same expression at all 
Mach numbers. 


The difference of the two remaining integrals gives the difference 
in the drag at Mq = 1 of a nonlifting wing and body having the same 
longitudinal distribution of cross-section area. Since the two integrals 
will not always cancel, the drag of the wing and body will, in general, 
be different. One can account in this way for the difference in drag of 
thin elliptic cone -cylinders and circular cane-cylinders disclosed pre- 


viously by integration of surface pressures. To show this, one must 
evaluate the integrals of equation (152) at the shoulder of the cone- 

cylinder (x = l) using the expressions for <p~ and <p 2 given in equa- 

W B 

tions (llO) and (93) . In the integration, the contour Cy extends on 
both sides of the y axis from -ml to 4 ml whereas the contour C B is 
a circle of radius 01 . Upon carrying out the indicated operations, one 
obtains the same result as that given previously in equation (122) in 


which the drag of the elliptic cone -cylinder is substantially less than 


that of the circular cone-cylinder. 


Special cases for which the drag of wing and body is the same .- 
Although it is important to note the difference in the drag of two bodies, 
it is perhaps even more important to know under what conditions the drag 
of the bodies is the same. If attention is confined to nonlifting cases 
so that qu is zero, the vortex drag vanishes , the condition for the 

equality of the drag of a wing and body having the same area distribution 
is that the contribution of the last two integrals of equation (153) cancel. 
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This condition is satisfied for certain large and important classes of 
shapes. One such class includes shapes that are cylindrical at the base 
since, then. 


and 




dn 


(15*0 


% = Db 


(155) 


Another includes many shapes that taper to a point at the rear since, 
then, both integrals again vanish. Other classes for which it is more 
difficult to specify the geometry include shapes for which the integrals 
have equal values different from zero. The latter case provides some 
interesting situations in which some members of a family of wings and 
bodies having the same longitudinal distribution of cross-section area 
have the same drag and others have a different drag. To be more specific, 
consider a low-aspect-ratio pointed wing having a straight trailing edge 
normal to the free-stream direction and smooth airfoil sections closing 
with a finite wedge angle at the rear (an example of such a wing is a 
tria n gular wing with biconvex profiles), and a body of revolution having 
the same area distribution as the wing. Application of equation (153) to 
this pair of bodies quickly leads to the conclusion that the drag of the 
body of revolution is Infinitely greater than that of the wing. This is 
apparent because the integral around the base of the wing is finite, 
whereas that around the base of the body is logarithmically infinite since 







Zn R 


L 


(156) 


and ds/ dx is finite and R is zero at x = Z . The infini te drag of 
this particular body of revolution is, of course, spurious' and is no 
doubt associated with the fact that the round stern is too blunt to treat 
with a theory of the slender-body type. On the other hand, there id no 
reason to believe that the pressure drags of the wing and body are the 
same. 


Since no corresponding difficulties occur at the base of the wing, 
let the drag of the above wing be compared with that of a second thin 
low-aspect-ratio wing having the same longitudinal area distribution. 
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For such a comparison, equation (153) must be replaced by the correspond- 
ing relation between the drag of two wings 



It is immediately clear that the two wings have the same drag if they have 
the same geometry at the base. This condition for the equality of drag 
of two bodies having the same area distribution has been arrived at pre- 
viously by Harder and KLunker (ref. l6) and by Beradt (ref. 29) by some- 
what dif ferent means. As is apparent from the preceding discussion, this 
condition is sufficient but not necessary. 

As a further example, consider the case where the geometry of the two 
wings is aff inely related, that is, for a constant chord, the second wing 
is derived from the first by simple stretchings of the y and z dimen- 
sions. For the present class of wings, having straight tra ilin g edges 
normal to the free-stream direction and airfoil sections closing with a 
finite wedge angle at the rear, each of the integrals of equation (157) 
can be written in the form 


Pz w1ta 


-n 


d(z/t) 


s 0 Td 


Then the product At of aspect ratio and thickness ratio is the same for 
both wings, although A and r individually may be different. It follows, 
moreover, from the fact that the ratio z/t of the wing ord ina tes to the 
maximum thickness is the same function of x/l and y/s Q for affinely 
related wings, that cp^ is the same function of y/s Q for both wings. 
Since s Q is proportional to A and At is the same for both wings, it 
follows that the two integrals of equation (157) have the same value, and 
both wings have the same drag. Inasmuch as it is only the conditions at 
the trailing edge that enter into the integrals of equation (157) > similar 
reasoning shows that the two wings also have the same drag if the condi- 
tion of affinely related geometry applies only to the cross-section shape 
and surface slopes in the x direction at the trailing edge. On the 
other hand, if the wings have merely the same longitudinal distribution 
of cross-section area, the simple relations just described between the 
various elements of equation (157) no longer hold, and the wings will, 
in general, have different drags. 
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Application to Uonplanar Problems 


Equation ( 90 ) expressing the relation between the perturbation 
potential for sonic flow about a thin low-aspect-ratio wing and that 
about a slender body of revolution has been derived on the assumption 
that the boundary conditions for the wing can be specified on a planar 
surface. The development outlined in sketch (f) suggests that the result 
can actually be extended to include more general classes of slender 
shapes. Accordingly, assume that equation (90) holds for cases involving 
non planar boundary conditions and let the results given in the prece ding 
sections for the drag at Mq = 1 of thin elliptic cones be extended to 
include slender elliptic cones of any eccentricity. The analysis proceeds 
identic ally to that for the thin elliptic cone, the only change being that 


<P ; 


2W 


must be recalculated. 


potential theory since 9 . 


% 


This is a simple problem in two-dimensional 
represents the potential associated with 


uniform growth of an ellipse, and the result, when evaluated on the cone 
surface, is 


% = 


U 0 tmx 

21 


to r“(i + -s_y 

.2 \ 2mls _ 


(159) 


Substitution of this result into equation (153) leads to the following 
relation for the drag 




m?t 2 Zn 


[i( 


1 + 2 ml) 


( 160 ) 


Comparison of these two expressions with the corresponding relations 
for thin cones given in equations ( 110 ) and ( 122 ) shows that they differ 

by the inclusion of an 



additional factor (l+t/ 2 mZ) 
in the more general result . 
Although the contribution 
of this term is of negli- 
gible importance for cones 
having t/m l small, it is 
vital for nearly circular 
cones, and indeed necessary 
to assure the equality of 
IU and Bg. when the ellip- 
tic cone Decomes a circular 
cone, that is, when t/mZ= 2 . 


2m . major axis 
vt = mror axis 


Sketch (p) 


point further, sketch (p) 
has been prepared, showing 
the variation with 2 ml/t, 
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or the ratio of major axis to minor axis, of the drag at Mq = 1 of two 
families of elliptic cones . All members of each family have the same 
longitudinal distribution of cross-section area. As indi cated, one 
family is defined by mt/Z = 0.02 and includes the circ ular cone-c ylind er 
having a semiapex angle 0 of 0.10 radian, and the other by mt/Z = 0.03 
and includes the circular cone -cylinder having 0 = 0.1225 radian. The 
solid line indicates the values computed using equation (160) , and the 
dotted line those computed using equation (122). In both cases, the 
drag Dg of the circular cone-cylinder is calculated from equation (127) . 
As would be anticipated, the solid and dotted lines coincide for thin 
cones, but they differ considerably for circular cones (2mZ/t = l) . More 
interesting, perhaps, is the extent to which the drag of a family of cone- 
cylinders having the same longitudinal distribution of cross-section area 
depends on the shape of the cross section. 

The procedures applied here to the elliptic- cone-cylinders can also 
be applied to many other cases, such as wing-body combinations, etc. For 
bodies having the same longitudinal distribution of cross-section area as 
a cone-cylinder one must merely determine the appropriate function for 
<p 2 jj and proceed in the same manner as for the elliptic cane-cylinders. 
For other bodies it is also necessary to have knowledge of either the 
theoretical solution or the experimental pressure distribution for sonic 
flow around a body of revolution having the same (or affinely related) 
longitudinal distribution of cross-section area as the given body. It 
should be remarked, however, that the extension to same of these problems 
involves the assumption that equation (90) applies to nonplanar cases. 


Comparison With Experimental Results 


In the remainder of this paper, experimental data will be presented 
and a comparison made with the predictions of sonic slender -body theory. 
Although these comparisons may not be ideal, since experimental data for 
Mq = 1 are only available for families of shapes that strain the assump- 
tions of the theory, they show remarkable agreement with the theory and 
help define the range for which the results may be expected to apply. 

Cone-cylinders .- The most informative class of bodies to investigate 
further with regard to comparison of theory and experiment is the cone- 
cylinder. This is because of the availability of not only the similarity 
rules, etc., but also the complete solution for the pressure distribution 
and flow field in the vicinity of such bodies. Experimental data are also 
ava i lable in reference 18 by Solomon for the pressure at several points 
on the surface of two rather blunt circular cone -cylinders at Mach numbers 
near unity. The tests were conducted on cone-cylinders having semiapex 
angles of 20° and 25° and at Mach numbers up to about O.96 . The corre- 
sponding pressures at Mo = 1 can be obtained only by extrapolation. The 
test Mach numbers are sufficiently high, however, that the local' Mach 
numbers on the body surface are virtually independent of the free-stream 
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Mach number . These pressures are plotted in sketch (q) together with the 
theoretical pressure distribution for slender cone -cylinders at Mq = 1 . 



The latter were computed using equation (112) together with the theoreti- 
cal pressure distribution for a circular cone-cylinder having semiapex 
angle Q = 0.10, see sketch (i) . With the exception of one point on the 
25 ° cone-cylinder, the theoretical and experimental values are in remark- 
able agreement, considering the bluntness of the canes. 

It would be very informative to make similar comparisons for cone- 
cylinders that are more slender, or that have noncircular cross sections, 
but the authors are unaware of any suitable experimental data. Studies 
involving bodies of revolution having area distributions that differ from 
that of cone-cylinders are handicapped at the present time by the lack of 
theoretical solutions for the transonic pressure distribution, and would 
have to be confined to the investigation of such items as the range of 
applicability of the similarity rules, the existence and lateral extent of 
the g(x) function, etc. 

Wings .- Since complete solutions for sonic flow around three- 
dimensional wings have not yet been obtained, the following discussion 
must be confined to cases in which experimental information is known for 
two or more wings or bodies having the same or affinely related longitu- 
dinal distributions of cross-section area. Probably the most extensive 
set of data of this type is that given in reference 20 for a large family 
of affinely related wings of rectangular plan form having MCA 63 AQXX 
sections. Since the results for Mq = 1 can be presented most concisely 
by using the variables suggested by the transonic similarity rules for 
wings of finite span (see, e.g., ref. 4), the experimental results for 
the zero-lift pressure drag and the lift-curve slope at Mq = 1 are 
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given in sketch (r) by plotting Gq/t 5 / 3 and ri/sc^ as functions of 
At 1 / 3 . As shown previously by McDevltt (ref. 20), these data confirm the 
statement provided by the similarity rules that the results so plotted 
should define a single curve for each aerodynamic quantity. 


DRAG AND LIFT AT M 0 = l 
RECTANGULAR WINGS- NACA 63A0XX PROFILES 



12 3 u I 

A(rr ACr-r 
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A666444432222L5II I I .5 
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Sketch (r) 


The curves representing the zero-lift pressure drag and the lift- 
curve slope have the same general form for high -aspect-ratio wings. The 
curves approach horizontal lines for the wings of larger aspect ratio, 
and the values for the lift and drag are not too different from the 
theoretical values given by Guderley and Yoshihara in references 30 and 31 
for two-dimensional sonic flow around double -wedge profiles. The curves 
approach straight lines through the origin for low-aspect-ratio wings and 

the experimental values for the lift-curve slope of wings having At 1 / 3 
less than about unity practically coincide with the theoretical values 
given by equation (1^3) . The corresponding theoretical values for the 
drag at Mq = 1 are not known. 

Same measure of the applicability of the theoretical results can still 
be derived, however, by examining the relation between the experimental 
drags of various wings. If the effects of the violation of the theoretical 
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requirements "that occurs at the round nose of the unswept leading edge of 
each of the present family of wings can he disregarded, the discussion 
following equation (158) applies and all low-aspect-ratio wings having a 
given longitudinal distribution of cross-section area have the same drag. 
Tnfl.RTmip.Vi as not many pairs of wings of the present family actually have 
the same area distribution, a more useful statement of the result is that 
the drag is a unique function of the area distribution. Since the area 
distribution of an affinely related family of wings can be specified by 
giving, for instance, the chord 2 and the ratio Sjn/l 2 (or its equiva- 
lent, the product of the aspect ratio and the thickness ratio) of the 
ma yiirnim cross-section area to the chord squared, it follows that the drag 
and geometry of the present family of wings are related according to 


D 

(Po U o 2 / 2 ^ 2 



t x (Ar) 


(161) 


where f x represents an unknown function of the indicated variables. 

This relation may be contrasted with that provided by the similarity rule 
that states 


Cp 

t 5/3 


=f 2 (AT 1/3 ) 


(p 0 U 0 2 /2)Sp- 


5/3 


(162) 


where Sp refers to the plan-form area. At first glance, the two 
relationships appear to bear only slight resemblance. It can be seen 
upon closer examination, however , that some of the apparent differences 
are superficial and of little or no significance. Thus, let equation (l6l) 
be rewritten as 


D = £oHo_ Z 2 f i ( At ) = ?E (At ) 2 f 3 (At ) 

or 

= AT 1/3 f 3 (Ar) ( 163 ) 

This appears to be the closest that the two results can be brought 
together without introducing additional restrictions or approximations. 

Both are now concerned with the same quantity, Cj)/r S//S , but equation (162) 
states that this quantity is equal to some unknown function of Ar 1/3 , 
whereas equation (163) states that it is equal to Ar 1/3 times some 
function of At. The only way in which both results can be universally 
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correct is for the functions f 2 and f 3 to be constants and not dependent 

on either At or Ar 1 ^ 3 . Both rules are not universally correct, however, 
since equation ( 163 ) is derived from transonic slender body theory and 
therefore can be expected to apply only to wings of small aspect ratio. 

From Idle foregoing considerations, one can conclude that the drag at 
Mq = 1 of the low-aspect-ratio wings of the present family must depend on 

the geometry in such a way that C^/t 5 / 3 is linearly proportional to A •W 3 . 
Examination of the drag data of sketch (r) shows that the experimental 

results exhibit precisely this trend for wings of At 1 ^ 3 less than about 
unity. An alternative interpretation of' this result is that the drag of 
a number of low-aspect-ratio wings of the present family all having the 
same chord, varies as the square of the frontal area; that is, 

D /(| p o u ° 212 ) deP 6115 - 8 on the square of Sj J~L Z - Further discussion of 

these and related points appears in reference 32 . 

It appears that the degree of correspondence between theory and 
experiment disclosed above for such extreme cases as rectangular wings 
of aspect ratios 3 and 4- must be attributed partially to the averaging 
influence of integration, and that the same close correspondence may not 
be found for more detailed quantities. For example, slender-body theory 
indicates that the lift on low-aspect-ratio rectangular wings at Mq = 1 
is concentrated along the leading edge. Although pressure-distribution 
measurements were not included in the test program reported in refer- 
ence 20, pitching-moment measurements were made from which the center-of- 
pressure position can easily be deduced. The results indicate that the 
center-of -pressure position at small angles of attack is within the first 
10-percent chord at Mq = 1 for each of the wings of aspect ratio l/2 or 1, 
but moves progressively rearward for wings of larger aspect ratio. Hence 

the range of Ar 1 ^ 3 for which theory and experiment agree may be expected 
to be less than that indicated by the integrated lift and drag results. 

On the other hand, application to wings of rectangular plan form imposes 
severe strain on the slender-body assumptions, and better agreement, or a 
wider range of applicability, might be anticipated with wings of other 
plan form, such as triangular. 

Wing-body combinations .- Several comparisons between the experimental 
zero-lift drags of wing-body combinations and bodies of revolution having 
the same longitudinal distribution of cross-section area were given by 
Whitcomb in reference 6 in connection with his discovery of the area rule. 
The bodies tested were of such geometry that the integrals of equa- 
tion ( 153 ) are zero and the drags of wing-body combinations and their 
equivalent bodies of revolution should be the same. The exper iment al 
results show excellent agreement in some cases, and lesser agreement in 
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other cases. These results will not he discussed further here since the 
experimental data are already analyzed from the point of view of equality 
of drag in reference 6. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Apr. 2, 1956 
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